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Q . Abstract 

I In any theory it is unnatural if the observed values of parameters lie very close to special 

values that determine the existence of complex structures necessary for observers. A natural- 
ness probability, P, is introduced to numerically evaluate the degree of unnaturalness. If P 
is very small in all known theories, corresponding to a high degree of fine-tuning, then there 
is an observer naturalness problem. In addition to the well-known case of the cosmological 
constant, we argue that nuclear stability and electroweak symmetry breaking represent sig- 
. nificant observer naturalness problems. The naturalness probability associated with nuclear 

stability depends on the theory of flavor, but for all known theories is conservatively esti- 
mated as Pnuc ^ (10~^ - 10"^), and for simple theories of electroweak symmetry breaking 
^ . -ffewSB ^ (10~^ - 10"^). This pattern of unnaturalness in three different arenas, cosmology, 

l/-^ I nuclear physics, and electroweak symmetry breaking, provides evidence for the multiverse, 

^!t| ■ since each problem may be easily solved by environmental selection. In the nuclear case the 

problem is largely solved even if the multiverse distribution for the relevant parameters is 
I relatively flat. With somewhat strongly varying distributions, it is possible to understand 

' both the close proximity to neutron stability and the values of mg and — m„ in terms of 

■ the electromagnetic mass difference between the proton and neutron, 6em — 1 ± 0.5 MeV. 

It is reasonable that multiverse distributions are strong functions of Lagrangian parameters, 
since they depend not only on the landscape of vacua, but also on the population mecha- 
nism, "integrating out" other parameters, and on a density of observers factor. In any theory 
with mass scale M that is the origin of electroweak symmetry breaking, strongly varying 
multiverse distributions typically lead either to a little hierarchy, v/M ~ (10~^ - 10~^), or 
to a large hierarchy, v <^ M. In certain multiverses, where electroweak symmetry breaking 
occurs only if M is below some critical value, we find that a little hierarchy develops with 
the value of /M"^ suppressed by an extra loop factor, as well as by the strength of the 
distribution. Since the correct theory of electroweak symmetry breaking is unknown, our 
estimate for Pewsb is theoretical. The LHC will lead to a much more robust determination 
of -PewsB; and, depending on which theory is indicated by the data, the observer natu- 
ralness problem of electroweak symmetry breaking may be removed or strengthened. For 
each of the three arenas, the discovery of a natural theory would eliminate the evidence for 
the multiverse; but in the absence of such a theory, the multiverse provides a provisional 
understanding of the data. 
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1 Introduction 

The Standard Model, with three gauge forces and three generations of quarks and leptons, has 
been laboriously constructed from decades of data involving collisions of particles at ever higher 
energies. It represents a triumph for symmetries, but leaves many questions about nature unan- 
swered. These questions fall into three groups. The first revolves around observations of the 
universe that are not described by the Standard Model, such as dark matter, dark energy, baryo- 
genesis and infiation. Secondly there is the question of electroweak symmetry breaking: while the 
Standard Model provides a mathematical description, it does not contain the physical descrip- 
tion. The quadratic divergence in the Higgs mass-squared parameter implies that electroweak 
symmetry breaking is determined by new physics at a scale M at which the Standard Model is 
embedded into a more fundamental theory. This extraordinary behavior is termed "unnatural" 
— the Higgs mass-squared parameter, m|, is very sensitive to changes in some parameter x of 
the fundamental theory, especially as M grows: 

dlnx m\ 

The final group of questions boils down to "Why the Standard Model?" These include why the 
forces arc the ones wc observe, why the quarks and leptons have the masses they do, and why 
there aren't other forces and particles. 

Since symmetries were the key to the Standard Model, it is natural to suppose that all three 
sets of questions will be answered by introducing further symmetries, and this has largely defined 
beyond the Standard Model physics for the last three decades. Yet there are three clouds on the 
horizon for this symmetry approach, one for each set of questions 

• There is no known symmetry to explain why the cosmological constant is either zero or of 
order the observed dark energy. It appears to take a very unnatural value. 

• For the Standard Model to be natural, M must be low; we should already have observed 
either signs of supersymmetry or signals of new interactions in precision electroweak observ- 
ables. 

• Some of the Standard Model parameters take values in a special region that happen to yield 
certain nuclear properties. 

It is possible that these clouds are the first indications of a breakdown in the power of sym- 
metry to determine fundamental physics. If our universe is just one among very many in an 
enormous multiverse, then observed universes will be those that contain certain complex struc- 
tures necessary for observers [1]. Such arguments from environmental selection can potentially 
solve the cosmological constant problem, and yield a statistical prediction for the dark energy 



2 



in observed universes [2]. In this paper we consider the extent to which nuclear stabihty and 
electroweak symmetry breaking provide evidence for environmental selection. 

Many physicists, however, are reluctant to countenance any form of anthropic argument. Why 
give up on the traditional, extraordinarily successful, methods of physics? How can we hope to 
understand the conditions for intelligent observers when we struggle even to define what life 
is? Why make the extraordinary leap of postulating an extra-horizon multiverse, which has the 
smell of a secular form of God? In short, many believe that appeals to the environment are an 
escape from true science and that, in the absence of data confirming the conventional symmetry 
approach, it would be better to change fields than to succumb to the philosophy of anthropics. 

The case of the cosmological constant demolishes these arguments. Traditional methods have 
not given any satisfactory understanding for why the cosmological constant is small. In contrast, 
the environmental argument not only explains why it must be small, but makes a statistical 
prediction for a non-zero value. Furthermore, this prediction requires essentially no understanding 
of what it takes to make an observer. Of course, the prediction does require a multiverse; but 
with a landscape of vacua provided by string theory [3] and generic inflation eternally creating 
universes [4] , a multiverse no longer appears theoretically unreasonable, and its consequences can 
be explored under the assumption that we are typical observers [5] . Indeed, since few doubt that 
environmental conditions explain why we find ourselves on the Earth as opposed to elsewhere 
in the solar system, the resistance to environmental reasoning is hard to understand. As for 
any physical theory, the real issue is whether sufficient evidence for it can be found to make it 
convincing. 

Even if there is a multiverse, and environmental selection accounts for the value of the cos- 
mological constant, it is far from clear what the implications are for the rest of physics. One 
extreme possibility is that all the other questions left open by the Standard Model will be solved 
within some symmetry framework, for example supersymmetric 5*0(10) with fiavor symmetries. 
In this case, the environmental solution of the cosmological constant problem simply allows us 
to ignore this special problem while solving everything else with symmetries — it justifies the 
model-building approach of the last 30 years. The other extreme is that the Standard Model is 
as far as we can get with symmetries, all of its parameters are either strongly selected by envi- 
ronmental effects or just reflect random typical values in the multiverse, and similarly for issues 
outside the Standard Model. Of course there are many intermediate cases between these two 
extremes. For example, there could be gauge coupling uniflcation within an SO{10) theory, but 
the flavor sector may involve sufficient scanning parameters that the electron, up and down quark 
masses are determined by environmental selection. The aim of this paper is to seek evidence 
for environmental selection both in the measured values of Standard Model parameters, and in 
electroweak symmetry breaking. 
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The cosmological constant illustrates an additional crucial point: environmental selection pre- 
dicts parameters that are, from the conventional viewpoint of symmetries, unnatural. This is a 
key point, because so much effort of the last 30 years has been expended in trying to understand 
how new symmetries could naturally lead to an unnatural effective low energy theory, namely 
the Standard Model. Perhaps the unnaturalness is a hint that symmetries are not the answer. It 
was realized some time ago that if only the mass parameter of the Standard Model scans, then 
environmental selection could solve the gauge hierarchy problem [6]. In the context of supersym- 
metry, differing assumptions about scanning parameters and their distributions could lead to a 
large hierarchy [7] or a small hierarchy [8]. These papers both illustrate that environmental se- 
lection leads to unnaturalness, and that measurements at the LHC could experimentally confirm 
the presence of unnaturalness in electroweak symmetry breaking. It is an open question whether 
the landscape favors a low or high supersymmetry breaking scale [9] . 

The first part of this paper, sections 2 to 5, are devoted to a consideration of the concept 
of naturalness, and its application to issues of nuclear stability. We are motivated by the belief 
that environmental selection may lead to precise predictions, and that it is important to evaluate 
the numerical significance of such predictions. We often hear it said that, in nuclear physics 
and elsewhere, nature exhibits "amazing coincidences" that give rise to life. How can these be 
precisely evaluated? In section 2 we formulate the concept of naturalness in a very general way, 
so that in section 3 we are able to identify naturalness problems associated with the existence 
of complex structures in the universe. In section 4 we study the stability of neutrons, deuterons 
and complex nuclei in the parameter space of the Standard Model, identifying how close our 
universe is to these stability boundaries. In particular, we find a closeness to the neutron stability 
boundary that has not previously been elucidated. The naturalness problems associated with the 
closeness to these stability boundaries are evaluated numerically in a variety of theories of flavor 
in section 5. We stress that the evidence for unnaturalness depends on the ensemble of theories 
being considered, but we find that there is an irreducible amount of unnaturalness no matter 
what the theory of flavor. 

In the second half of the paper, we investigate the consequences that follow if this unnat- 
uralness arises from environmental selection in a multiverse. It is here that the real utility of 
our new formulation of naturalness becomes apparent. In section 6 we argue that problems of 
unnaturalness are indeed solved by environmental selection, and show how the amount of unnat- 
uralness is connected to the probability distribution of the multiverse. We argue that evidence 
for environmental selection will increase if symmetry arguments fail to solve an accumulation of 
naturalness problems. In the case of nuclear stabilities, we show in section 7 how this leads to 
multiverse predictions for the masses of the electron, up and down quarks. In section 8 we embed 
the Standard Model Higgs sector in a generic theory of electroweak symmetry breaking at mass 



4 



scale M. Allowing parameters of this sector to scan, as well as parameters of the Standard Model, 
we consider environmental selection from nuclear stability boundaries. We find that sharply vary- 
ing multiverse distribution functions generically lead to both large and little hierarchies between 
M and the scale v of electroweak symmetry breaking. From the viewpoint of the multiverse, a 
discovery at the LHC that electroweak symmetry breaking requires some amount of fine-tuning 
would not be surprising. A similar analysis on electroweak symmetry breaking is performed in 
section 9, assuming that the relevant boundaries are the phase boundary of electroweak symme- 
try breaking rather than the nuclear stability boundaries. We again find that unnaturalness in 
electroweak symmetry breaking is expected for sharply varying multiverse distribution functions. 
In section 10 we study connections between environmental selection for the cosmological constant 
and for electroweak symmetry breaking. For example, we argue that the driving force on the 
multiverse for unnaturalness in electroweak symmetry breaking could arise from the distribution 
for the cosmological constant, the connection being through weakly interacting massive particle 
(WIMP) dark matter. Finally, our conclusions are given in section 11. 

2 New Definition of Naturalness Problems 

The concept of naturalness has often been the driving force for finding fundamental mechanisms 
in nature. Arguments for naturalness are often phrased in terms of the sensitivity of low energy 
Standard Model parameters q to variations of the parameters in the more fundamental theory 
ttj. A simple measure of naturalness is then given by A = \d\n.Ci/ d\naj\, with a large value of A 
signaling a lack of naturalness [10]. This definition, however, could miss some of the important 
aspects of the naturalness problem. Here we introduce a new, quantitative definition of natural- 
ness that can be applied in much more general situations. In particular, this allows us to identify 
certain classes of naturalness problems that have not been quantitatively defined. Our definition 
is also free from some of the problems existing in the simplest definition of naturalness based on 
the sensitivity of parameters. 

First of all, it is very important to notice that in talking about naturalness, we are dealing, 
either explicitly or implicitly, with an ensemble in which parameters of the theory are varied 
according to some definite distribution. Consider, for example, that the Higgs mass-squared 
parameter, is given by the difference of two mass scales of order the fundamental scale, 
m\ — Ml — M|. We say that the theory is unnatural if |m^| ^ Ml.M^. This statement, 
however, already assumes that it is unhkely for Mi and M2 to be very close or for both to be very 
small; more specifically, the distribution of possible values for Mi and M2 is assumed to be almost 
structureless in the M1-M2 (or Mf-M|) plane. This illustrates that the concept of naturalness is 
closely related to the distribution of parameters in an ensemble. 
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A particular member of an ensemble is unnatural if it has parameters very close to special 
values that are not explained by the symmetries of the theory. The parameters are special if 
some physical property arises that is not a generic feature of the members in the ensemble, or if 
they separate regions of parameter space that have differing generic features. Depending on the 
property considered, we encounter various classes of naturalness problems, whose solutions could 
point to various different mechanisms in nature. With this definition, the degree of unnaturalness 
is given by how close the parameters arc to the special values, which gcncrically form a special 
hypersurface in multi-dimensional parameter space. As we will see below, we can quantify this 
degree in terms of the distribution of parameters within the ensemble. In this section, after 
presenting a new definition of naturalness, we apply it to well-known situations. In the next 
section we use it to introduce new types of naturalness problems. 

2.1 Definition 

Let us start describing our precise definition of naturalness by introducing the distribution func- 
tion f{x) for an ensemble. The function is defined such that the number of members with 
parameters Xi {i = 1, - ■ ■ ,N) taking a value between Xi and Xi + dxi is given by 



Here, Xi represent continuous parameters of the theory, e.g. masses and coupling constants.^ The 
overall normalization of / becomes relevant if we are interested in the total number of members, 
which will be the case when we discuss relative likelihoods of different theories. For the present 
purpose of discussing naturalness of a given theory, however, the overall normalization of / is not 
important. 

The distribution function depends on the choice of Xi. What variables should we choose as Xj? 
In general we can choose any variables as Xi, depending on the context. To discuss naturalness 
of the low energy theory, however, it is often most convenient to take "fundamental" parameters, 
such as masses and coupling constants of the ultraviolet theory, as Xj. Now, suppose we have only 
one such variable x. If the observed value of x, Xo, is very close to a special value x, the level of 
unnaturalness is given by the following "naturalness probability" 



"'^In principle, we can define a distribution function that also has discrete labels representing, e.g., particle 
content and symmetries of the theories. Here we consider a distribution function for each theory that has a 
definite symmetry, matter content, number of spacetime dimensions and so on, and restrict its arguments Xi to be 
continuous parameters associated with that particular theory. 



dN" — f{xi, X2, - ■ ■ , xn) dxidx2 ■ ■ ■ dxN- 



(2) 



P 



lTf{x)dx 



(3) 
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where x^ai-n < 2; < Xmax gives the range of x values in the theory under consideration.^ Here, we 
have included only members on one side of the special point, but depending on the situation, a 
factor of 2 should be added to the numerator to include members on both sides of the special 
point. With this definition, P <^ 1 signals the existence of a naturalness problem. 

The form of the distribution function is modified if we redefine the variable x. For example, 
if f{x) = l/x for a variable x, the change of the variable x' = \nx can make the distribution 
function flat for x': f{x') = 1. It is often convenient to go to the basis in which the distribution 
function is constant. In that basis, the naturalness probability of Eq. (3) becomes 



Xn 



X 



X-n 



Xj] 



(4) 



which is simply given by the distance between Xo and x divided by the available parameter space.^ 
As should be the case, this quantity becomes smaller as Xo approaches x. 

The naturalness probability of Eqs. (3, 4) can be extended to the case of multiple parameters 
Xi. The precise deflnition depends on the dimensionality of Xi and the special surface consisting 
of Xi- If the co-dimension of the special surface is 1, we can choose a; to be a linear combination 
of Xi perpendicular to the surface, and then use the deflnition of Eq. (4) (in the basis where the 
distribution function is constant). In the case that the co-dimension is higher, we must use an 
appropriate generalization of Eqs. (3, 4) deflned using multi-dimensional volumes, rather than 
simple distances. A useful deflnition, in the basis where f{xi) is constant, is 



(5) 



where a — 1, ■ ■ ■ ,n runs over variables whose observed values, Xa,o, are close to the special values, 
Xa, and c„ = n'^^'^/T{n/2 + 1) is the volume of the unit ball in n dimensions. Here, the expression 
in the numerator, Vn, is a measure of the accidentally small volume of parameter space required 
for Xa,o to be close to Xa, while that in the denominator, Vn, is the total volume of parameter 
space available. The restriction of the variables to those with \ {xa,o~^a)/{xa,max — Xa,m.m)\ <^ 1, 
e.g. < 1/5, is important to avoid obtaining P <^ 1 simply as a result of high dimensionality of 
the parameter space. In the case that the special surface in question arises as an intersection 
of co-dimension 1 surfaces, we may restrict the volume of the numerator to be one side of the 
co-dimension 1 surfaces, depending on the situation. This definition reduces to that of Eq. (4) in 
the special case oi n — 1. 



2 With f{x) = for range of the integration in the denominator can be taken from 

—00 to +00. Practically, if f{x) has a sharp drop-off, it is useful to restrict the range of the variable by a sharp 
cutoff, as Xmin and Xmax iu Eq. (3). 

^This interpretation is also possible for a general distribution function f{x) if we consider f{x) to be a sort of 
metric, or volume factor, in parameter space. 
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Figure 1: Illustrations of the definition of the naturalness probability P in multi-dimensional 
parameter space. 

The definition of Eq. (5) is illustrated for a 2-dimensional parameter space (xi, X2) in Fig. 1(a). 
The curve C represents special values of these parameters corresponding to some physical phe- 
nomenon, and the curve C similarly represents special values for some other phenomenon. If the 
observed values (xi^o, 0:2,0), denoted by the dot, are close to the intersection point (xi,a;2), the 
naturalness probability is given by Eq. (5), with a = 1,2. The numerator represents the area of 
the shaded region, and the denominator represents a much larger area corresponding to the range 
of the parameters. Here, we have restricted the area of the numerator to one side (the lower side) 
of the curves C and C", anticipating an application in later sections. 

In Fig. 1(b), it is apparent that the observed point is much closer to C than C"; but the 
definition of Eq. (5) does not capture the additional unnaturalness associated with this. (The 
area of the numerator in this expression is represented by the dotted line.) In this situation, 
we must first consider the naturalness probability Pc associated with C, using the definition of 
Eq. (4), where the axis x is taken normal to C and passing through the observed point at a; = Xo- 
We then take the coordinate x' along the curve C and consider the naturalness probability Pc' 
associated with C in this coordinate. The resulting naturalness probability is then P = PcPc, 
which is given by the area of the shaded region in the figure divided by that corresponding to the 
range of the parameters. This tends to zero as the observed point approaches the curve C, as 
expected. In general, if we find several unnatural features that have a hierarchy in their degrees of 
unnaturalness, we can obtain the correct estimate for the naturalness probability P by considering 
it as a product of several naturalness probabilities, each associated with an unnaturalness of some 
fixed degree. The decomposition of P can be made along the lines presented here. 
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Figure 2: The contour of m\/Ml = ±1, ±0.1, ±0.01, ■ ■ ■ in the M1-M2 plane. The special line of 
m| = is visible at Mi = M2. The observed value of M2 — Mi + O(10~^^) is denoted by the 
little dot for an arbitrary value of Mi = 1.2M*. 

2.2 Illustrations 

Let us now illustrate the use of our definition in the case of the conventional gauge hierarchy 
problem. Consider that the Standard Model is embedded into the fundamental theory at a high 
scale of Ml ^ M2 ~ and that the Higgs mass-squared parameter is given in terms of Mi 
and M2 by m\ = Ml — M|. In general, we expect that the distribution function is roughly fiat 
in terms of Mj {i = 1,2) and that the range of parameters is < Mj < 0(M*).*^ In Fig. 2, 
we plot the contour of m\ in the M1-M2 plane in units of M*. We clearly see that there is a 
special line in this plane, Mi = M2, where = 0. The physical phenomenon making this 
line special is that members very close to this line have a low energy effective field theory that 
contains a scalar excitation whose mass is hierarchically smaller than M*. For ^ 10^^ GeV, 
the observed value of m| ~ —0(100 GeV)^ is denoted by the httle dot for an arbitrary value of 
Ml = 1.2M*. We find that it is located very close to the special line compared with the expected 
range of Mj. In fact, this is always the case regardless of the value of Mi we choose, and is a 
manifestation of the lack of naturalness in the Standard Model. The problem is that the special 
line with m| = is not special in terms of the symmetry structure of the theory. The points with 
= are no more symmetric than those with 7^ 0. By choosing a variable x to be a linear 
combination of Mi and M2 perpendicular to the special line, and then using the definition of 
Eq. (4), we obtain P 10~^^. A similar analysis can also be made for the cosmological constant 

*Here, we have assumed that Mi are mass parameters associated with "fermions" (including mass parameters 

in supersymmetric theories). The argument, however, is not affected if they are associated with scalars, in which 
case the distribution is expected to be roughly flat in Mf with the range — O(M^) < < 0{M^). 
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problem. In general, the cosmological constant A is given by the sum of terms of order taking 
the form of MfM'j. The observed value of A is then very close to the special hypersurface with 
A = 0, and the points with A = are no more symmetric than those with A 7^ 0. The physical 
property making this hypersurface special is that members very close to this hypersurface allow 
an arbitrary large observable universe. From Eq. (4), the degree of unnaturalness we obtain in 
this case \s P ^ 

Here we note that our definition docs not suffer from the problem existing in the simplest 
definition of naturalness based on the logarithmic derivative, A = \d\nci/ dhiaj\. Suppose that 
a dimensionful parameter is given by two dimensionless constants and b as ^ = e~^'^^^^^^, 
where is the fundamental scale, and that a natural range for and b in the fundamental theory 
is g'^b — 0(0.1). We then naturally obtain /x/M^, = e~*-'(^°°°); there is nothing unnatural with this. 
The simplest definition, however, gives A = |(91n/i/(91n = 0(1000), signaling (incorrectly) the 
existence of unnaturalness. Our definition does not lead to such a fake signal, since any point 
with g% — 0(0.1) is generic — there is no special value of g% that can be singled out as x in 
Eq. (3, 4). 

We now see how theories beyond the Standard Model solve the gauge hierarchy problem in 
our language. As we have seen, the Standard Model embedded into more fundamental theory 
at M* 10^^ GeV leads to the probabihty of having a small weak scale v fa 100 GeV to be 
P {v/M^Y 10~^^. This is because the Higgs mass-squared parameter, m|, is given by 
the sum of various contributions of order M^, with the distribution function expected to be 
roughly fiat in Mi (or M?). The distribution function, however, varies from one theory to another. 
For example, if the weak scale arises from dimensional transmutation associated with some gauge 
coupling g becoming strong, as in technicolor [11] and certain supersymmetric theories [12], then 
V ~ exp(— 87r^/|6|(?^). Here, b is the one-loop beta function coefficient for the new gauge force 
and g^ = g{M^). The distribution function is expected to be roughly fiat in the variable g1,^ 
giving an approximately fiat distribution for l/ln(f/M*). To illustrate this situation, we plot in 
Fig. 3 the contour of ml = Mf - M| with Mi = cxp(-87rV|&|£/*) and M2/M1 = r in the 
\b\g'^-r plane. This can be regarded as a simplified model for the situation in, e.g., dynamical 
supersymmetry breaking. Typically, there is considerable uncertainty in the range of \b\gl. Here 
we take < \b\gl < 0(10). As is seen from the figure, the observed value of v can correspond to 
a completely generic (not special) point in the plane, which is indicated by the little dot for an 
arbitrary value of r = 1.4. This, therefore, provides a solution to the gauge hierarchy problem; 
see Fig. 2 for comparison. The important point here is that different theories can lead to very 
different distribution functions f{x) for physical parameters Xi. A distribution flat in v"^ is the 
origin of the gauge hierarchy problem, while one flat in In v (or 1/ In v) provides the solution. The 

^The conclusion below does not change if the distribution is flat in g*, and so on. The distribution, 

however, could affect arguments on certain naturalness problems; see section 4. 
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Figure 3: The contour of m\/Ml — 10 ^" with n = 3, 4, 5, ■ ■ ■ in the |6|5'^-r plane. The observed 
value of f ~ y'— m| ^ 100 GeV can correspond to a completely generic point in the plane, which 
is indicated by the little dot for an arbitrary value of r = 1.4. 

situation with the cosmological constant is similar, although the problem is much more severe 
because all known theories give distributions that are essentially fiat in A, while none are known 
that are fiat in In A (or 1/ln A). 

2.3 Generality of the framework 

So far, we have illustrated our definition of naturalness by looking at naturalness problems that 
have a conventional description: the observed values of parameters are (very) close to special 
values where the sensitivity of a low energy parameter to high energy ones diverges, for example 
(91nm^/(91nM(, (91n A/SlnMj oo. This is one class of naturalness problems, which we may 
call "fine-tuning naturalness problems." The real power of our formalism, however, lies in the 
fact that we can discuss many different types of naturalness problems in a unified manner, simply 
by extending what we mean by "special values" for parameters. Any physical property that is 
not a generic feature of the members in the ensemble is a candidate for identifying special values. 
There can be many classes of naturalness problems, depending on the property considered. The 
closeness of the observed values to the special values, signaled by P ^ 1, can then be used as 
evidence for a new mechanism to understand such accidents. 

One application of this general idea is to use naturalness arguments as evidence for the pres- 
ence of some symmetry beyond the Standard Model. For this purpose it is often convenient to 
take parameters of the low energy theory to be Xj. Imagine that the observed values Xi^o of N di- 
mensionless Standard Model parameters Xi {i — 1, • • ■ , N) take values close to the special surface 
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that defines a symmetry relation 

S{xi) = 0. (6) 

If the symmetry were in fact absent, nature would be described by some member of an ensemble 
giving the Xi parameters distributed according to f{xi), which is not generically concentrated on 
the surface of Eq. (6). After redefining the parameters to make the distribution function flat, 
we can introduce an axis x normal to the symmetry surface that passes through the observed 
point at Xi — Xi^o- The degree of unnaturalness is then given by Eq. (4). The closer the observed 
point is to the special value x on the symmetry surface, the less probable it becomes that nature 
is described by this ensemble, i.e. by the theory without the symmetry. A small value for P 
can thus be taken as evidence that the symmetry is indeed present in nature, at least in an 
approximate form.^ It could be that Xo — x is dominated by the experimental uncertainty. In this 
case improved experiments may provide further evidence for the symmetry. A well known example 
is provided by a grand unified symmetry, such as S'f/(5), that gives a symmetry relation between 
the Standard Model gauge coupling constants S{gl, g^-, g1) = [13]. If in fact there is no unified 
symmetry, the observed values of a, and sin^^^^ give a naturalness probability of P ~ 0.1 
in an ensemble without supersymmetry, and P ~ 0.01 in an ensemble with supersymmetry [14]. 
Here, a, and 6w are the fine structure constant, the effective QCD coupling, and the Weinberg 
angle, respectively. This lack of naturalness can then be regarded as evidence for some form of a 
unified symmetry in theories with supersymmetry. 

In the next section we introduce classes of naturalness problems that arise from the existence 
of complex structures. We will see that the formalism developed here elucidates the identification 
of these naturalness problems. As in other naturalness problems, successful solutions to these 
problems could lead us to find new mechanisms or dynamics in nature, which we consider in later 
sections. 

3 Complexity and Observer Naturalness Problems 

The definition of naturalness introduced in the previous section allows us to identify new classes 
of naturalness problems. A member in the ensemble is unnatural if it has parameters unusually 
close to "special" values; but clearly there are many reasons that parameters could be "special." 
We frequently stress special values that lead to a large hierarchy of mass scales, but in this section 
we consider special values that lead to the existence of relatively long-lived complex structures, 
such as nuclei, stars and galaxies. 

^Here we assume that conversions from low energy parameters to high energy ones, e.g. through renormalization 
group evolution, do not significantly affect the size of P. This is generically a good assumption. 
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3.1 Complexity naturalness problem 

As the parameters Xi vary, moving from one member of an ensemble to another, suppose that 
a physical threshold is crossed that is crucial for the existence of some complex structure. This 
defines a special surface in the parameter space 

C{x,) = 1, (7) 

that divides the volume of parameter space into two regions, one that supports the complex 
structure, C{xi) < 1, and one that does not, C{xi) > 1. In general this surface is not one 
with enhanced symmetry. Therefore, a member in the ensemble having parameters unusually 
close to this surface has a "complexity naturalness problem." The degree of unnaturalness is 
given numerically by Eq. (3), where the single variable x is normal to the "complexity surface" 
of Eq. (7), where it takes the value x, and the particular unnatural member in question has a 
nearby value for this parameter, Xo- 

One caveat is that we have assumed that the physical threshold relevant for the existence of 
the complex structure is sharp in the parameter space. This can be verified in any particular case, 
and is certainly true in the examples discussed in the next section. For example, the stability of 
any particular nucleus gives a sharp boundary corresponding to values of the coupling strengths 
and quark masses that lead to a surface of zero binding energy. In general a lack of sharpness 
is due to the time evolution of complex structures in universes corresponding to the different 
members of an ensemble. As parameters vary from one member to another, complex structures 
could gradually become less stable. For example, the stability of large scale structure is not 
completely sharp — as the cosmological constant is gradually increased only the regions with 
larger statistical fluctuations in the density perturbations are able to collapse. Still, the relevant 
parameter space for the cosmological constant spans over 100 orders of magnitude, and over this 
space the transition for the existence of large scale structure is very sharp. 

Another caveat is the hidden assumption that these complexity surfaces are not distributed 
so densely throughout the entire parameter space that a typical member in the ensemble is 
expected to be close to one or more surfaces. There may indeed be many complexity surfaces; 
for example, in the Standard Model there are several hundred relatively stable nuclei, each with 
its own complexity surface in an ensemble that contains the Standard Model. However, if the 
relevant parameters — the Yukawa couplings, y, the weak scale, v, and the QCD scale, Aqcd — 
vary by many orders of magnitudes in the ensemble, then these complexity surfaces will all be 
tightly clustered in a "complexity zone," yv ~ ckAqcd ~ O(0.01)Aqcd- A complexity naturalness 
problem now arises because this zone is itself small compared with the entire volume of parameter 
space. Most members in the ensemble lie in voids far from the complexity zone, and the closeness 
of the observed parameters to one or more of the complexity surfaces implies that the member 
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describing our universe lies in a special region. There may be several complexity zones for nuclear 
physics; for example ones with four, five or six quark flavors lighter than the QCD scale, but 
in certain parameter directions each zone will be small. In the case of a single variable, the 
degree of unnaturalness can be taken to be P = Axz/ (xmax — a^min) where Ax^ is the width of the 
complexity zone. If the density of complexity surfaces is very high it might be that a member 
in the complexity zone is very close to some surface; but this could be a reflection of the density 
of the surfaces rather than any additional unnaturalness beyond that of being in the complexity 
zone. 

In summary, a theory possesses a complexity naturalness problem if, in the parameter space 
of an ensemble, the member describing our universe lies very close to a surface corresponding to a 
physical threshold that allows the existence of some relatively long-lived complex structure. For 
the problem to exist, there should not be many surfaces of a similar character distributed densely 
and almost uniformly over the parameter space; otherwise, the closeness to one of these surfaces 
would simply be a generic phenomenon for the members in the ensemble. 

3.2 Observer naturalness problem 

Now we introduce another, closely related, naturalness problem. Complex structures are required 
for the existence of "observers." This implies that some complexity boundaries may also act 
as boundaries that divide the parameter space of the ensemble into those members that may 
support certain observers and those that can not. These boundaries are harder to define than 
general complexity boundaries, since we are unable to give a precise definition of an observer. 
Moreover, the parameter region may not simply be divided into the two regions "with" and 
"without" observers — the expectation value for the number of observers may in general be 
a complicated function over the parameter space, with the value significantly varying across a 
complexity boundary. Nevertheless, since the changes of the expectation value across some of 
these boundaries are expected to be drastic, caused by drastic changes of complex structures, for 
certain purposes we may approximate this function to be step-like. This leads to the concept of 
the "observer boundary" 

0{xi) = 1, (8) 

which divides the parameter space into one "with" observers, 0{xi) < 1, and one "without," 
0{xi) > 1. In general, there are many complexity boundaries, Ca{x) — 1, that are relevant for 
the existence of observers. The region of parameter space that allows observers, O, is then given 
by the common set of Ca{x) < 1 for all a, and the observer boundary, 0{xi) — 1, is the border of 
this region. As more boundaries are added as relevant ones for observers, the region O shrinks, 
but the number of observer boundaries docs not generically increase. This is the crucial difference 
between complexity and observer boundaries. For illustration of these boundaries in an example 
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(b) 



Figure 4: Characteristic situations for complexity and observer boundaries in 2- dimensional pa- 
rameter space. Complexity boundaries that are and are not relevant for the existence of observers 
are depicted by solid and dashed lines, respectively. The shaded region indicates an observer re- 
gion O, which (a) may or (b) may not be a small region around the observed point, which is 
denoted by the dot. The observer boundary, 0{xi) = 1, is given by the border of O and is 
represented by the thick solid line. 

of 2-dimensional space, see Fig. 4. 

While the identification of observer boundaries suffers from some ambiguities, there are certain 
advantages to focusing on these boundaries rather than on general complexity boundaries. In 
general there can be several disconnected regions, Oj, that can support observers. However, since 
the existence of observers undoubtedly requires certain complex structures, we can be convinced 
relatively easily that these regions are not distributed densely throughout the entire parameter 
space. For example, if we take the existence of stable complex nuclei to be one of the required 
conditions for observers, then we find that almost the entire region with v/Aqcb ^ 10"^, with the 
other parameters of the Standard Model fixed, is outside Oj [6], regardless of the existence of 
any other physical thresholds. Increasing the size of the up or down Yukawa coupling more than 
a factor of a few, with the other parameters fixed, also leads to a sterile universe in which both 
proton and neutron cannot be stable inside nuclei. The situation is similar for other conditions. 
For example, the requirement for the existence of complex structures in the universe, such as 
galaxies, makes the entire region with A/Q^T^^ ^ 1 outside Oj [2], where Q is the primordial 
scale of density perturbations and Teq the temperature of matter radiation equality. These imply 
that the observer regions Oj fill only a small fraction of the parameter space. They are small 
islands in the entire parameter space of the ensemble, and so the observer boundaries are not 
distributed densely and uniformly over the parameter space. 
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The sparseness of Oj allows us to focus on the region O (c {Oj}) that contains our observed 
point. The unusual closeness of the observed parameters to the relevant observer boundary (the 
boundary of O) then signals unnaturalness, raising an "observer naturalness problem." The degree 
of unnaturalness can be quantified using the definition given in the previous section. The estimate 
for the degree will be conservative if we use only crucial requirements for the existence of observers, 
such as the existence of nuclei, stars, and galaxies, since this will give a large allowed region for 
observers, O. Imposing more and more subtle conditions, such as the existence of carbon, oxygen 
and so on, will only decrease O and therefore increase the degree of unnaturalness. A relevant 
question is whether we find an observer naturalness problem in the Standard Model and beyond, 
keeping only relatively robust conditions. This question will be addressed in the next section. 

4 Evidence for Observer Naturalness Problems 

In this section we study if there exist observer naturalness problems in the Standard Model and 
beyond. In general it is easier to study the existence of an observer naturalness problem than that 
of a complexity naturalness problem, since it is easier there to see that the relevant boundaries 
are not uniformly and densely distributed over the parameter space. In order to identify and 
quantify the observer naturalness problem reliably, we take the following approach. We first list 
complexity boundaries that are candidates for composing the observer boundary, in the relevant 
parameter space of the Standard Models of particle physics and cosmology. These boundaries 
have effects of various degrees on the environment, and wc take only the ones giving fairly drastic 
effects. This leads to a large observer volume O and hence to a conservative estimate for the 
degree of an observer naturalness problem. The actual degree of unnaturalness could only be 
more severe if other complexity boundaries are important for the development of observers, since 
this would reduce O. 

4.1 The relevant parameters 

Let us consider some ensemble that leads at low energies to the Standard Model of particle physics 
(taken to include neutrino masses) and to the standard cosmological model. Different members of 
this ensemble give different values for the parameters of the Standard Model, x^u, and cosmology, 
Xcos- Some members may give values for this set x = {xsm^^Jcos} that are so distant from the 
observed values, Xo, that the corresponding physics and astrophysics is completely different from 
that observed — for example, if five quark fiavors were lighter than the QCD scale, or if the baryon 
asymmetry were of order unity. Here we restrict our discussion to a subset of the ensemble in 
which variations in the parameters about Xo arc limited. In fact, we are interested in determining 
whether small deviations in x from Xo could dramatically change certain relatively stable complex 
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structures (atoms, stars ...) that we observe. 

The physics and astrophysics of these complex structures depends on only a subset of x; for 
example, variations in the bottom quark mass, or its mixing to the charm quark, by a factor of 
2 have no effect on the complex structures of interest. Hence, we now restrict x to include only 

xsM ■■ a, as, Vu, Vd, yt, Ve, ml, (9) 

and 

Xcos ■ Mpi, Teq, Q, rjB, A. (10) 

Here, yu,d,t,e ai'e the Yukawa couplings for the up, down and top quarks and the electron, the 
Higgs quartic coupling, m\ the Higgs mass-squared parameter, Mp\ the reduced Planck mass, 
and rjB the baryon asymmetry. In this 13-dimcnsional parameter space there are special surfaces 
that represent complexity boundaries 

CAix) = 1, (11) 

such that the complex structure A differs drastically from one side of the boundary to the other. 
Several comments are in order for our choice of the parameter set x. 

• \h and m\ can be traded for the clcctrowcak vacuum expectation value, v, and the Higgs bo- 
son mass. We include \h and yt because they typically play an important role in clcctrowcak 
symmetry breaking. For example, in the Standard Model there is a boundary corresponding 
to the existence of a electroweak symmetry breaking vacuum with v -C Mpi [15] 

C,{Xh,yt)^l. (12) 

• The SU{2) gauge couphng of the Standard Model is omitted from Eq. (9), since charged 
current weak interactions at low energies are described by v. Neutral current, and therefore 
the weak mixing angle, play little role on relevant complex structures. 

• We have omitted the strange quark Yukawa coupling yg for simplicity. It does play a role in 
nuclear physics. 

• The only parameter from the lepton sector is ye- Other lepton sector parameters could 
affect TjB via leptogenesis [16]. Since the source of rjB is unknown, we prefer to list it as an 
independent cosmological parameter. 

• The 13-dimensional parameter set could be further reduced, since only certain combinations 
appear in Ca of Eq. (11). For example, only dimensionless combinations appear, so that Mpi 
could be removed by using it as the unit of mass. For the complexity boundaries arising from 
atomic and nuclear structures, v appears only in the combinations m/ = yfV {f — u,d,e) as 
fermion masses, so that for these boundaries xsm can be reduced to 

, mu rud rUe , . 

-'iQCD Aqcd Aqcd 
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where Aqcd is the QCD scale, which takes a value of 100 MeV in our universe. For 
naturalness arguments, however, we often prefer to use the more basic set of Eq. (9). 

• In a theory with two Higgs doublets, such as the minimal supersymmetric standard model, 
the basic set of Eq. (9) should be expanded. For most purposes, we simply have to replace the 
electroweak vacuum expectation value v by two vacuum expectation values; for type-II two 
Higgs doublet theories, we have Vu for the up-type Higgs doublet and va for the down-type 
Higgs doublet. 

4.2 The relevant complex structures 

A crucial question for the complexity boundaries is what are the relevant complex structures. Our 
interest in a particular boundary depends on how important the corresponding complex structure 
is for explaining the structure of the physical world. Consider two extremes. The size of the 
cosmological constant, relative to the matter density of the universe when density perturbations 
become non-linear, leads to a boundary that determines whether large scale structure forms: 

A 



A member on one side of the boundary will lead to formation of galaxies, while one on the other 
side leads to an inflating universe of isolated elementary particles. This clearly has important 
effects on the basic structure of the universe. As a second extreme example, consider the complex 
mesons with a h quark constituent. As the h quark mass, mb, is increased above the W boson mass, 
the h quark decays so rapidly that B mesons cease to exist. This boundary of mi,/my\/ ^ 1, 
however, has a negligible effect on our environment, which is why we did not include the 6-quark 
Yukawa coupling, i/b, in xsm- 

Somewhat arbitrarily, we divide the relevant complexity boundaries into three classes accord- 
ing to how dramatic the environmental change is across the boundary: 

• Catastrophic boundaries change our universe into one that is essentially unrecognizable. In 
addition to Eq. (14), and Eq. (12), we would include the case that electroweak symmetry is 
broken dominantly by the Higgs potential 



V 



.1. (15) 

^^QCD 

This is required so that the baryon asymmetry of the universe is not washed out by the 
sphaleron effects [17]. We also consider that the absence of any complex nuclei is catas- 
trophic. In the simplified case that the only parameter that is varied from its Standard 
Model value is v, this boundary is [6] 

^^TT ~1- (16) 

IO^Aqcd ^ ' 
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For V larger than this boundary, the only stable nucleus is either p — uud or A++ = uuu. 

• Violent boundaries separate members where a crucial complex structure of our universe is 
absent. For example, across the boundary 



where nip and m„ are the masses for the proton and neutron, the neutron becomes stable 
and hydrogen unstable. Such a neutron-stable world would not have dense astrophysical 
objects fueled by nuclear energy release, such as main-sequence stars in our universe [18]. 
Another example of violent boundaries is that of vanishing deuteron binding energy 



Across this boundary, no nuclei form during big bang nucleosynthesis, so that the universe 
protonizes (or neutronizes). Stars could only burn via exotic triple proton reactions, with 
extremely high central densities. 

• Substantial boundaries separate members where a crucial complex structure of our universe 
is drastically changed. For example, across certain boundaries stars may exist but are very 
different from those we see. For example, if the pp reaction pp De^v is not available to 
ignite stars, then protostars would collapse to a higher temperature before igniting via the 
pep reaction, ppe^ Dv. A more substantial change to stars occurs if the deuteron is beta 
unstable, D ppe^u; stars could still burn, but only by using the helium produced during 
big bang nucleosynthesis. The existence of a stable diproton, the ^He nucleus, would also 
change stellar nuclear reactions, shortening lifetimes of hydrogen burning stars. Changes of 
some nuclear energy levels, controlled by quark masses and a, could also lead to substantial 
changes of the abundances of various nuclear species, such as carbon and oxygen. 

Each boundary gives a surface of special values for the parameters, x, which does not corre- 
spond to a surface of an enhanced symmetry. Among the ones listed, we find that most boundaries 
are clustered around a zone yu,d,eV ~ ctAqco ~ O(0.01)Aqcd, where our universe also resides. 
This strongly suggests the existence of a complexity naturalness problem. Here, however, we focus 
more on the observer naturalness problem, which arises if the member describing our universe lies 
unusually close to the observer boundary. There are certain ambiguities in identifying which of 
the complexity boundaries compose the observer boundary. An important point, however, is that 
by using only the boundaries that certainly have disastrous effects on the environment, we can 
be on the conservative side in evaluating the existence of an observer naturalness problem. For 
this reason, we take only boundaries that are catastrophic or violent, rather than just substantial, 
to compose our observer boundary. We also select only the boundaries that allow us to derive 
reasonably accurate values for x, allowing a reliable estimate of the naturalness probability P. 
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4.3 Stability boundaries for neutrons, deuterons and complex nuclei 

Following work by others, we focus on the boundaries across which we lose complex nuclei, 
Eq. (16), neutron instability, Eq. (17), and the deuteron, Eq. (18). We take yd, v/Aqcd 
and a to be our parameters Xj. The choice is motivated by the expectation that the relation of 
these parameters to those in the ultraviolet theory is relatively direct. 

We now represent the boundaries of Eq. (16, 17, 18) in terms of the deviations of Xi from 
the observed values Xi^o- Let us first discuss the neutron stabihty boundary of Eq. (17). The 
instability of a neutron, or equivalently the stability of hydrogen, requires 

m„ — rup — nie > 0, (19) 

where we have neglected the neutrino mass. The neutron-proton mass difference, m„ — nip, arises 
from both the strong isospin violating effect, Sd^u, and the electromagnetic contribution to the 
proton mass, 5em' — nip — S^-u ~ <^em- In our universe, S^-u — 2.26 ± 0.51 MeV [19] and 
(^EM = 5d-u — {f^n — f^p) ~ (2.26 ± 0.51) — 1.29 MeV. To a first approximation, these quantities 
scale as 5d-u oa nid — niu and 5em oc a, so that 

nin-nipc^ ^(2.26 MeV) (0.97 MeV). (20) 

'md,o - niu,o OLo 

Here we have taken the central value for bd-u, and the variables with and without the subscript 
o represent, respectively, the values in our observed universe and those in an arbitrary member 
of the ensemble. In terms of the parameters xi, Eq. (19) can be written as 

4- { y^-y^ - 0.23^') j^P^ > 0.56^. (21) 

Here, we have neglected a small dependence of Aqcd on v as well as logarithmic evolution of the 
Yukawa couphngs. For two Higgs doublet theories,^ the Yukawa couplings in Eq. (21) should be 
replaced as 

Vu Vu sin p, yd yd cos p, y^ ye cos P, (22) 

yu,o yu,o sin Po, yd,o yd,o cos Po, ye,o ye,o COS po, (23) 

where v = \Jv'^-\-v^ and tan/3 = v^jvd- 

We next discuss the boundary of Eq. (18). This boundary corresponds to the stabihty of a 
deuteron under strong interactions, which requires 

Bd — mp + nin — niD < 0, (24) 



''Here and below we assume type-II two Higgs doublet models when we discuss two Higgs doublet theories. 
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where m/j is the deuteron mass. The binding energy Bd depends on the quark masses as 

B.-B.,„ = -af^^^4^-l). (25) 

\mu,o + md,o J 

where Bd,o — 2.2 MeV is the observed deuteron binding energy. The parameter a is uncertain, 
but it is estimated in Ref. [6] to be a ~ (1.3 - 5.5) MeV using models of nucleon binding. This 
aUows us to write Eq. (24) in terms of 

Vu + yd v/AdCD ^ ^ ^ 2.2 MeV ^^g) 



yu,o + yd,o {v/A^qcd)o 

For two Higgs doublet theories, the Yukawa couplings must be replaced as in Eqs. (22, 23). 

We finally consider the condition for the existence of complex nuclei. The stabihty of complex 
nuclei requires the energy release for the (3 decay, n — > pe~P, to be smaller than the binding 
energy of nuclei per nucleon, Eki^: 

mn-mp-me^ -E'bin, (27) 

where we have neglected the neutrino mass. Precisely speaking, E'^in varies with a nucleus, 
receiving contributions both from nuclear forces and the Coulomb repulsion between protons. 
To a first approximation, however, E'bin can be regarded as the same for all nuclei and being 
controlled purely by Aqcd, taking the value of E^n — 8 MeV in our universe.^ Substituting this 
into Eq. (27), and using Eq. (20), we obtain 

' y'-y- 0.23^1//^ ^4.6 + 0.56^. (28) 



0.77 \yd,o ~ yu,o ' ye,oJ {v/Aqcd)o ' ' a 
The expression should be modified according to Eqs. (22, 23) for two Higgs doublet theories. 



4.4 Observer naturalness problem in the Standard Model and beyond 

The region inside our observer boundary, O, is defined by Eqs. (21, 26, 28), with modification 
by Eqs. (22, 23) for two Higgs doublet theories. To visuahze the overall shape of this region, in 
Fig. 5 we depict the three boundaries of Eqs. (21, 26, 28) in m„-md-me space. The plots in the 
left column adopt a logarithmic scale in m„ and m^, while those in the right a linear scale in m„ 
and TTid- The fine structure constant is fixed to be a = ao, and we have used the leading-order 
chiral perturbation value of mu,o/nid.o = 0.56 to draw these plots. (The qualitative features of the 
plots are not affected if we vary S^-u and mu,o/nT'd,o in the range of 6ci-u — 2.26 ± 0.51 MeV [19] 
and mu,o/'<^d,o = 0.3 - 0.6 [20].) From the figure, we find that the observed parameters are close 

^We neglect the dependence of i^bin on m„ and rud, as this will not affect our conclusions. 
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mu/mu,i 



iTiu/mu, 




Figure 5: The location of the observer boundary in mu-md-rrie space. The neutron and complex 
nuclei boundaries of Eqs. (21) and (28) are depicted by solid lines (below and above, respectively). 
Dashed lines represent the deuteron boundary of Eq. (26) (for a = 5.5, 2.2 and 1.3 MeV from 
below). The observed point is represented by little dots. 
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to the three boundaries, especially when viewed on a logarithmic scale. In particular, they are 
very close to the neutron stability bound of Eq. (19). This implies that we have an observer 
naturalness problem. 

The degree of unnaturalness P depends on a theory, since the relation between {m„, m^, me-, a} 
and the fundamental parameters, as well as the form of the distribution function, depend on the 
ensemble we consider. However, we can still make a conservative estimate on P, based on the 
observation that the naturalness probability P obtained by fixing all the fundamental parameters 
except for one generally satisfies P ^ P. In the present context, we can fix y^, Vd, Ve, ol and 
^QCD; as well as tan/3 for two Higgs doublet theories. The conditions of Eqs. (21, 26, 28) then 
lead to 

0.5 ^ - < 2, (29) 

where we have used a — 2.2 MeV for illustrative purpose. The probability of v falling in this 
range, P^,, then gives a conservative estimate of the naturalness probability P. Alternatively, we 
can fix V instead of Aqcd- In this case, we have 

0.5 <, -^^^ ^ 2. (30) 

The probability of Aqcd falling in this range, P^^^-^, can also give a conservative estimate of P, 
since P ^ min{P^, Paq^^^}. Which of P„ and Paqcd use is determined by which is easier to 
estimate and, in case both are easy to estimate, which gives a stronger bound on P. 

Let us first consider an ensemble in which the distribution of v is flat in v^, as in the Standard 
Model embedded into more fundamental theory at a high scale. In this case, we obtain an 
extremely small value of P^ ^ v^/M^, where M* » Vg is the cutoff scale. For M* 10^^ GeV, 
we obtain Py 10"^^. The origin of this small value, however, is precisely the existence of the 
gauge hierarchy problem. We thus see the existence of a severe observer naturalness problem in 
these theories, but it is hard to disentangle from the flne-tuning naturalness problem. 

We therefore focus on theories in which the conventional gauge hierarchy problem is solved. 
In these theories, the distribution of v is expected to be flat in Inw or l/\nv within an ensemble. 
To estimate P^, however, we need to know the range of v, which depends on how v is related to 
the fundamental parameters of the theory. To avoid model dependence coming from this, we here 
consider Paqcd; instead of P^. The value of the QCD scale, Aqcd, is determined from the strong 
gauge coupling constant at the cutoff scale, (73 (M*), through renormalization group evolution. 
The probability Paqcd can then be estimated if (i) the theory between Aqcd and M^,, together 
with the vahic of M^,, is specified, and (ii) the distribution function for g^{M^), including the 
range of g^{M^), is given. ^ Specifying the theory between Aqcd and is important because 

^Here we neglect higher order effects, which are expected to be small. 
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Figure 6: The value of gl{M^) as a function of Aqcd/Aqcd,o for non-supersymmetric theories 
(left) and supersymmetric theories (right). The range of Eq. (30) is depicted by the vertical lines, 
while the corresponding range of gl{M^) by the horizontal lines. 

the existence of colored states whose masses are associated with v may increase the correlation 
between v and Aqcd, enhancing -Paqcd- The distribution and the range of g^{M^) are also 
important. For example, if the distribution of g^i^M^) were flat in l/g^{M^) within the range 
l/g^{M^) ^ b^/lGn^, the QCD scale is distributed almost flat in IuAqcd for all Aqcd ~ M^. 
Here, 63 is the one-loop beta function coefficient for g^ evaluated at M*. In this case, we would 
obtain -Paqcd estimated using Eq. (30) to be inflnitely small. To avoid this unphysical conclusion, 
we can introduce an arbitrary cutoff c on the distribution, l/gl{M^,) ^ c. Alternatively, we can 
consider that the distribution of g3{M^) is flat in g'^{M^). In this case, restricting the range to be 
gl{M^) ^ c' gives a finite answer to Paqco- The value of -Paqcd depends on c', but we can make 
a reasonable conjecture on the value of c', e.g. c' 1 or Idn'^/bs. 

In Fig. 6, we plot the value of gl{M^) as a function of Aqcd/Aqcd,o obtained using one- 
loop renormalization group equations for non-supersymmetric theories (left) and supersymmetric 
theories (right). In non-supersymmetric theories, is taken to be 10^^ GeV, which is the scale 
where the SU{3) and C/(l) gauge couphngs almost meet. (This is the scale where the three gauge 
couphngs would almost meet if five Higgs doublet fields or a vector-like fermion with the Higgs 
quantum numbers were introduced at the weak scale.) In supersymmetric theories, we take to 
be the unification scale, 2 x 10^^ GeV, and the scale of superparticle masses to be msusY — 1 TeV. 
Dependence of the results on these parameters, however, is weak. From the figure, we find the 
probability of Aqcd falling in the range of Eq. (30) to be 
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for non-supersymmetric theories and 

Paq,„ ~ 0.038 ^ ^ 26-^ i (32) 

for supersymmetric theories. Here, we have taken the distribution function to be flat in g^{M^). 
The value of c' is uncertain, but we expect c' ^ 1. The numbers in Eqs. (31, 32) depend on a in 
Eq. (25), through the dependence of the lower value of Eq. (30) on a. For a — 1.3 MeV (5.5 MeV), 
0.016 and 0.038 in these equations become 0.020 and 0.046 (0.012 and 0.029), respectively. Note 
that the analysis here provides the most conservative estimate for the level of an observer natu- 
ralness problem existing in theories beyond the Standard Model (in which the conventional gauge 
hierarchy problem is solved, so that Py can be larger than Paqcd)- 

-P^-Paqcd- (33) 

In theories where the gauge hierarchy problem is not solved, the degree of the observer naturalness 
problem is much more severe, P ^ P^, -C 1. We thus find that the bounds on P derived here 
provide evidence for an observer naturalness problem in the Standard Model and beyond. 

To obtain a conservative estimate on P, we have used here only P ^ mm{Py, Paqcd}- Physi- 
cally, the case with P ^ min{Py, P\^^^} corresponds to the situation where we have the "best" 
model of flavor, i.e. the observed quark and lepton masses in units of the weak scale are repro- 
duced essentially without any free parameter. This is, however, not the case for many theories of 
flavor, which will have P < min{Py, Paqcd}) ^ discussed in the next section. 



5 Naturalness Probabilities and Theories of Flavor 

The observed values of the first generation masses, mu,d,e-, are special: Fig. 5 shows that they are 
close to the boundaries of neutron, deuteron and complex nuclei stability. To judge whether this 
is likely to be accidental we need to evaluate the naturalness probability, but this depends on the 
theory of flavor. We begin by considering ensembles of theories in which a and i'/Aqcd are flxed, 
but the flavor observables xp have some probability distribution f{xp). 

In the Standard Model the relevant flavor observables are the Yukawa couplings: xp — yu,d,e- 
While they are all very small, of order 10~^ - 10~^ at unifled scales, they are all close to the 
maximal values allowed by the neutron, deuteron and complex nuclei boundaries. Of particular 
important, the observed values are especially close to the neutron (hydrogen) stability boundary. 
While Uf, is about a factor of 3 from this boundary, in the yu-Ud plane the distance from the 
observed point to this boundary corresponds to a variation in the coupling only of 
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where z = mu/md and A = m„+mp+me)/5d_u evaluated in our universe. Here, the central 

value is obtained with z = 0.56 and Sd-u = 2.26 MeV. The electromagnetic interaction raises 
the proton mass above the neutron mass by about an MeV. Given that quarks can have masses 
up to 100 GeV or more, it is remarkable that the rud — niu mass difference just overcompensates 
the electromagnetic shift to make the hydrogen mass only ~ 0.78 MeV smaller than the neutron 
mass. The deuteron and complex nuclei stability boundaries are also close by, corresponding to 
changes in the Yukawa couplings by factors of ^ (1.5 - 3) and ^ 2, respectively. 

It is important to numerically evaluate these accidents; for illustration, we consider two simple 
distribution functions. If f{yu,d,e) ai'c fiat and non-zero in the range of yu,d,e from to 1, then 
the corresponding naturalness probability for this ensemble, using Eq. (5), is Pp ~ 10^^^; while 
if /(logj^o yu,d,e) are fiat and non-zero in the range of log^o yu,d,e from —6 to 0, then Pp ~ 10~^. As 
expected from Fig. 5, the volume of parameter space closer to the special points on the observer 
boundary than the measured point is very small compared with the expected total volume of 
parameter space in the ensemble. Prom the viewpoint of neutron instability and deuteron and 
complex nuclei stability, the Standard Model description of flavor is highly unnatural. 

Theories of flavor that go beyond the simple Yukawa coupling parameterization of the Standard 
Model typically involve further symmetries, such as flavor or unified gauge symmetries. While 
the set of flavor observables xp'va these theories can be smaller than that in the Standard Model, 
all known models do involve free parameters. If a theory could be found in which yu,d,e are 
precisely predicted, i.e. if variations in yu,d,e are less than the distances to the special points on 
the observer boundary when xp vary, then such a theory would have Pp ~ 1. However, if yu,d,e 
vary significantly (and independently) a,s xp vary in the ensemble, as in most theories of flavor, 
then the theory is likely to have Pp (much) smaller than of order unity. The best hope for a 
significant improvement in naturalness is then to obtain a successful prediction for the ratios of 
the first generation Yukawa couplings — symmetries that successfully predict the yd^ and y^. 
ratios could show that our universe lying so close to the observer boundary is just accidental. 
This is, however, not so easy, especially because of the extreme closeness of the observed point to 
the neutron stability boundary. Moreover, since the observer boundary involves the combinations 
yu,d,eV /-^QCB, the normalization of yu,d,e is coupled to the value of f /Aqcd- Below, we do not take 
v/Aqcb as fixed, but consider it to vary in the ensemble as discussed in the previous section. 

In certain theories of flavor, for example with Abelian flavor symmetries, the smallness of y/ 
(/ = u, d, e) follows from a single small parameter of the theory 

yf = CfeP^, (35) 

where e <^ 1 while c/ are of order unity. The powers pf are not free parameters of the theory, but 
result from a certain judicious choice of Abelian charges for all the quark and lepton fields [21] . The 
ratios of the Yukawa couplings can then be just numbers of order unity for a fixed value of e (and 
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tan (5 for two Higgs doublet theories). Let us consider, for example, that the coefficients c/ vary in 
the range < c/ < 2c/,o with a flat distribution on a linear scale, or in the range log^o Cf.o ^ 1/2 < 
logio — logio ^f,o'^^l'^ with a flat distribution on a logarithmic scale. In this case, the extreme 
closeness of the neutron stability boundary gives Pp 1/18 and 1/36, respectively, for a fixed 
value of e (and tan/?). In addition to this, there is the issue of the normalization of the masses 
mu,d,e, arising from a variation of v/Aqcb- Assuming that the distribution function is flat in 
g^[M^), this leads to a further reduction of the probability at least by a factor of P ^ 1/35 
and 1/15 for non-supersymmetric and supersymmetric theories, respectively. (These numbers are 
obtained by requiring that w/Aqcd should be within a factor of ^ {2 - 3) from the observed value 
so that its variation can be absorbed into changes of Cf, without a large extra cost in Pp. The 
range of gl{M^,) here has been chosen to be < g1{M^) < 1.) Possible variations of e (and tan/3) 
may or may not lead to further reductions of the probability. We thus conclude that theories 
with Abelian flavor symmetries considered here have 

P^PfP^ (10-3-10-2), (36) 

from physics of flavor alone. 

To increase P from flavor physics, it is necessary to predict ratios of yu,d,e to better than 
the factor of 2 we used in the above estimates. A more elaborate charge assignment for the 
quark and lepton fields under an Abelian flavor symmetry is not sufficient, since it will leave 
the relative coefficients of order unity undetermined. A non- Abelian flavor symmetry, or unified 
gauge symmetry, could improve P by successfully predicting yu '■ Vd '■ Ve as Clebsch factors. For 
example, : j/^ : ^/e = 1 : 2 : 1 at the unified scale agrees with inferred values of mu^d,e, providing 
the eigenvalues are not affected much by mixing to heavier generations. A simple theory of fiavor 
incorporating such a relation could lessen the significance of the observer naturalness problem of 
Eq. (36). However, uncertainties from QCD arc too large for us to know if the relation really puts 
us sufficiently close to the neutron stability boundary, so that the problem is simply an accidental 
consequence of symmetry. Moreover, in conventional theories of flavor we expect simple Clebsch 
factors to apply to the heaviest generation not the hghtest. This is because the hghtest generation 
typically gets contributions from nontrivial matrix diagonalizations, so that the relation would 
have to involve more than one generation. For example, the inferred values of mu,d,e are consistent 
with unified boundary conditions y„ = ye and yd — 0, provided the down quark mass arises from 
diagonalizing a symmetric 2x2 matrix that fully accounts for the Cabibbo angle. This requires, 
however, experimental inputs of the strange quark mass and the Cabibbo angle, so that it does not 
provide a real symmetry solution to the problem. We find it rather difficult to have a symmetry 
understanding for yu '■ yd'- ye, and even if we had, unnaturalness still arises from the coincidence 
between v and Aqcd, given by Eqs. (31, 32). We conclude that, in the absence of a convincing 
and significant progress in flavor physics, we have unnaturalness associated with flavor at the 
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(b) (c) 

Figure 7: The situations that lead to an observer naturalness problem. The observed point is 
denoted by the little dot, while the region outside O is shaded. 

level given by Eq. (36). 

6 Environmental Selection 

We have seen that theories of particle physics and cosmology are likely to have observer naturalness 
problems, signaled by P <^ 1. The problems arise in an ensemble if the observed point is close to 
the observer boundary. The situations in which this happens can be classified into the following 
three cases: the observed point is (a) a typical point within a small observer region O, (b) close 
to the boundary of a small observer region O, and (c) close to the boundary of a large observer 
region O, which may or may not be a closed region in the parameter space. These cases are 
depicted schematically in Fig. 7. 

The heart of the observer naturalness problem is the coincidence: why is the observed point 
so close to the observer boundary, which is a very special region in the parameter space? To solve 
this problem, we must make the observer boundary really special — the existence and location of 
the observer boundary should somehow affect the process of selecting a member in the ensemble 
as the one describing our universe. This leads us to consider environmental selection, effects 
sensitive to the existence of an observer. In this section we study how environmental selection on 
a multiverse can solve observer naturalness problems in general, and under what circumstances 
(and in what sense) we can identify evidence for it. We also discuss possible implications of this 
solution in identifying the correct theory describing our universe. 

The nature of the observer naturalness problem is dramatically altered if the members in 
an ensemble are physically realized in a multiverse. Until now, the ensemble has been a useful 
fictitious mathematical device to study problems of naturalness; now we assume that each member 




28 



of the ensemble is physically realized as a universe in spacetime. The fundamental theory of 
nature is assumed to have a huge number of vacua, the landscape of vacua, so that physics at 
low energy may be described by many possible effective theories. In fact, we are interested in 
a subset of these theories, T, that lead to physics described by the parameters Xi discussed in 
section 4, with the parameters Xi taking different values in different universes. We assume that 
it is possible to define a distribution function of the multiverse, /^(xj), such that, averaging over 
the entire multiverse, we obtain (xj) = J fxixi) Xidxi. We stress that this is the average value 
on the entire multiverse, independent of whether universes contain observers. If the population 
mechanism, including any relevant volume factors, were independent of Xi, then frixi) would 
represent the distribution of vacua in the theory underlying the landscape. However, in general 
it includes both relevant population and volume factors. ^° We first focus on a single theory, so 
that we can drop the index T until later subsections. 

Environmental selection on the multiverse involves two key points that give a central role to 
measurements made by observers 

• Universes that have parameters outside the region O do not contain observers, hence values 
of Xi outside O cannot be measured. When discussing naturalness of the observed universe, 
we should not be asking questions about the entire multiverse but only about the observer 
region O. 

• The number of observers in universes in O with ,Xj in the region Xj to Xi + dxi is given 
by f{xi) n{xi) dxi = f{xi) dxi. Each universe is to be weighted by an observer distribution 
n{xi), the factor associated with the number of observers that develop in a universe with 
parameters Xi. We make no attempt to define an "observer." 

These two points are really two aspects of a single selection process, with n{xi) defined over 
the entire space of Xj. We, however, find it useful to consider that n(xj) rapidly drops to zero 
at certain observer boundaries, so that the region outside these boundaries does not affect the 
naturalness of a multiverse. 

In the approximation that the observer boundary is sharp, the naturalness probability of 
Eq. (3) is then replaced by 

/I7(^) dx 



- r ./ N , , (37) 

where the integral in the denominator implies that x is integrated only in the region O. The 
expression for higher dimensional x space is also obtained similarly. We stress that here and 
throughout the rest of the paper 

f{xi)^~f{xi)n{xi). (38) 
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For subtleties in defining the distribution function in the multiverse, see e.g. [22] and references therein. 
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The entire region of parameter space of the multiverse is now irrelevant; the only question is 
whether we are typical observers in universes in O. To calculate / it is necessary to know both 
the landscape of vacua and the population mechanism. On the other hand, n is independent of 
the landscape and, with sufficient understanding of the physical environment for observers, could 
be calculated in principle from the low energy theory. We take the practical viewpoint that both 
/ and n are unknown, and hence give ourselves the freedom to assume any reasonable smooth 
distribution for /. It is then clear that distributions can be found that make Po ~ 1, solving the 
observer naturalness problems whether of the form of (a), (b) or (c) of Fig. 7. A more detailed 
description of each of these three types of solution is given in section 6.1. In section 6.3 we 
consider the relative probability of different theories, T, solving an observer naturalness problem. 

It is quite clear that there are many origins for a significant Xi dependence of f{xi) — 
f\xi) n{xi) in the observer region O. For / these are the distribution of vacua in the under- 
lying theory as well as relevant volume and population factors, while for n a variation in Xi could 
change the density of galaxies, stars, nuclei and so on that are relevant for observers. In sec- 
tion 6.2 we investigate a further effect: it could be that / or the shape of O depends on more 
variables, x^, than the ones we focus on, Xa, and that when these extra variables are integrated 
out the resulting effective distribution acquires further Xa dependence. Finally, a crucial issue 
is how to evaluate evidence for environmental selection. In section 6.4 we argue that, for any 
observer naturalness problem, a low value for the naturalness probability, P, provides evidence 
for environmental selection. More precisely, for the known set of (simple) symmetry theories, T, 
the evidence is determined by the largest value of Pt 

P^^ = maxr{Pr}. (39) 
If Pmax ~ 1 there is no evidence, but as Pmax decreases so the evidence becomes more substantial. 



6.1 Three manifestations of environmental selection 

Using a set of Lagrangian parameters Xi of some low energy theory, an observer naturalness 
problem might appear in the three ways illustrated in Fig. 7. In each case environmental selection 
may be at work, but the description of the solution to the problem is different. 

Suppose we find an observer naturalness problem of the form of (a) in Fig. 7. This suggests 
that the distribution function, /, is (approximately) constant in the Lagrangian basis Xj, since the 
replacement of P — > P© then completely eliminates the observer naturalness problem. Assuming 
a constant /, the naturalness probabihty Pq of Eq. (37) becomes 



Po 



Ax 



(40) 
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P <C 1 Po ~ 1 

Figure 8: Replacement of the probability P ^ Pq, illustrated for the case of (a) of Fig. 7. 

where Ax is the range of x in the observer region O. Similarly, for multiple parameters the 
naturalness probability of Eq. (5) is replaced by 



P 



O 



(41) 



where is the volume of parameter space in the observer region O. The elimination of 

observer naturalness problems by the replacement P Pq is illustrated in Fig. 8. The enormous 
increase in the probability here arises simply because all of the universes that lie outside the 
region O are cut out of the denominator. 

The situation in Fig. 8 is the case of the overall picture for the cosmological constant A. For 
conventional naturalness, consider an ensemble with a distribution function / that is constant in 
A, which is the case in most theories. The probability that a member of this ensemble has a small 
value of Ao near the special value zero is -P(A) = Aq/M^, where M^, is the fundamental scale. On 
the other hand, environmental selection, arising from a multiverse with / and n assumed fiat in 
A, replaces A^^ax ~ by AA ^ Q^T^^, giving 

P.(A) « (42) 

The vast majority of universes have a huge cosmological constant; but that is irrelevant because 
they contain only a dilute gas of elementary particles undergoing inflation. We should cut all 
these universes out of the measure, and consider only those where large scale structure forms, 
allowing the possibility of complex observers. This eliminates (or greatly ameliorates) the problem 
associated with the cosmological constant: -P(A) ^ 1 — Pc'(A) ~ 1. A more refined analysis 
includes the effects of nontrivial n{A) near the observer boundary [23]. 
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missed boundary 




Po <C1 ? 



Po 



~ 1 



Figure 9: Illustration of the situation in which a missed boundary (shown dashed) leads to a 
fictitious problem of Po ^ 1 after environmental selection. 

Let us now consider observer naturalness problems in the forms of (b) or (c) of Fig. 7. In 
these cases, even after selection we are apparently left with some residual naturalness problem 
Po <C 1 in some particular theory T. What does this imply? Does this mean that the observer 
naturalness problem in these forms cannot be solved by environmental selection? 

The answer is clearly no. To understand what can be going on, it is important to realize that 
our knowledge can often be very incomplete. Imagine that we have assumed some distribution 
function / for an ensemble and found an observer naturalness problem in the form of (b) or (c) 
of Fig. 7 in the basis where / is constant. This seems to imply that the observer naturalness 
problem cannot be solved by a simple cutout procedure illustrated in Fig. 8. However, do we 
really know that we have identified all the relevant complexity boundaries, and thus identified the 
correct observer boundary? It could be that we have missed some relevant boundary and, after 
taking it into account, Po may become order unity. This is illustrated in Fig. 9 for the case of 



Another possibility is that our initial assumption on the distribution function may not be 
correct. In practice, when we consider a theory, we start by assuming a "natural" distribution 
function /(xj), which is often taken to be constant in "Lagrangian parameters" Xj. However, there 
are many reasons why the distribution function in this parameter basis may not be constant even 
approximately. First, the Lagrangian parameters Xi of the low energy theory may be functions 
of the Lagrangian parameters yj in the fundamental theory, = Xi{yj), and the distribution 
function may be constant in yj rather than Xj. Second, in the multiverse picture, a mismatch 
between the naive and true distribution functions may also result from the population mechanism 
depending on Xj, or may be induced by environmental selection via a nontrivial n{xi). Finally, 



(b) of Fig. 7. 
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Figure 10: Illustration of the situation in which naive and true distribution functions differ 
significantly. In the basis where the true distribution function is constant (the mathematical 
basis), the observed point is a typical point in O, and the parameters chosen, yi, are in general 
nontrivial functions of "Lagrangian parameters" Xi. In the basis spanned by Xi (the physical 
basis), the distribution function is nontrivial, giving a "probability force" field, defined in Eq. (43). 

environmental selection may depend not only on Xi but also on other variables, inducing an extra 
dependence of / on Xj, as we will discuss in the next subsection. Hence, the distribution function 
f{xi) may have a strong dependence on the Lagrangian parameters Xi of the low energy theory, 
and in this case the situations of (b) or (c) in Fig. 7 actually correspond to typical observers with 
Po ~ 1, expected from environmental selection. 

The situation in which naive and true distribution functions differ significantly is illustrated 
in Fig. 10 for the case of (c) of Fig. 7. (For the discussion here we assume that all relevant 
parts of the observer boundary have been correctly identified.) There are two (equivalent) ways 
to describe this situation. One is to go to the basis in which the true distribution function is 
constant. In this basis, which we may call the "mathematical basis," the observed point is a 
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typical point in O, and the parameters chosen, y^, can in general be nontrivial functions of the 
"Lagrangian parameters" Xi. In practice, the transformation from Xi to Ui should be inferred from 
the observed data and a calculated observer boundary. Another way to describe the situation 
is to keep using Lagrangian parameters, or parameters that have the most direct or intuitive 
physical meaning at low energies, Xj. In this basis, which we may call the "physical basis," we 
have a nontrivial distribution function /(xj). The distribution function is peaked towards the 
region outside (9, so that a typical observer in O measures x,i close to the observer boundary. 

In the physical basis, Xj, it is very convenient to introduce a "probability force," which is 
simply the gradient of the probability distribution: 



where V = {d/dxi,d/dx2,- ■ ■)■ Within the region O, this vector field gives an indication of 
what values of physical parameters are most likely to be measured in the multiverse. If the 
vector field indicates a fiow towards an observer boundary, observers should be hving close to 
the corresponding boundary. An example of the probability force field is depicted in the physical 
basis plot of Fig. 10. 

6.2 Effective distributions from "integrating out" parameters 

What are the possible origins of the probability force? One is a nontrivial prior distribution 
function /prior, which arises from a mismatch between the physical and mathematical bases, as 
discussed previously. In practice, however, this is not the only source of the probabihty force. The 
key is that when we discuss physical predictions in the multiverse picture, we typically choose 
a subset, Xa, of the entire parameter set, Xi, and study probabilistic predictions in the space of 
these chosen parameters Xa- An important point here is that when we focus on Xa, we cannot 
simply ignore the other variables Xf, {i — {a,b}); they must be integrated out. This provides an 
effect on the distribution function defined in the subspace Xa, and thus modifies the probability 
force in Xa space. 

To illustrate this idea, let us consider the simplest example of a constant /prior, i-e. the case 
in which the physical and mathematical bases coincide. In this case the distribution of observers 
is given by 



where c = /prior is a constant, and the probability force in Xi space is zero, F = 0. This, however, 
does not mean that the probability force is zero if we consider only a subset of the variables, e.g. 
3^1, • • • , xn^i. When we focus on the variables Xi, - • • , Xjv-i (as we focused only on y„, yd, y^ 
v/Aqcd and a out of the 13 parameters of Eqs. (9, 10) in sections 4 and 5), we must integrate 
out the other variable xn- Now, in the multiverse picture with environmental selection, the 



F = V/, 



(43) 



dj\f — cdxi ■ ■ ■ dxN, 



(44) 
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only relevant universes are those within the observer region O, implying that integrals should 
be performed only over O. The domain oi Xn integration is then determined by the observer 
boundary, which is generically written as • • • , x^v-i) < xn < x^^{xi, • ■ ■ , xjv-i)/^ giving 

dj\f = fesixi, xn-i) dxf- dxN-i, (45) 

where fcs{xi, ■ ■ ■ ,xi\f_i) = c{xY^^^{xi, ■ ■ ■ ,xiy_i) — a;™™(a;i, ■ ■ ■ , X7v_i)}. The probability force 
defined in xi - ■ --xn^i space is then obtained by Eq. (43) with / replaced by fcs- — ^feS, 
where V = {d/dxi, ■ ■ ■ , d/dx^-i)- This generically gives a nontrivial force, Fcg 7^ 0. 

In general, when we consider physics of the landscape in the subspace of (a = 1, ■ ■ ■ , n), we 
must integrate out the other variables Xi, (6 = n + 1, ■ ■ ■ , A^) to obtain the effective distribution 
function, and thus the probability force, defined in Xa space. The domain of integration is given 
by the observer region C, leading to 

/eff(a;i, • • • , Xn) = j '" /prior(a:i, • • • , Xn) dXn+1 ■ ■ ■ dxN, (46) 

where /prior = /priori IS the prior distribution function defined in the entire parameter space Xi, 
whose Xb dependence could arise from both /prior and n. The important point here is that the 
effective distribution function fes in Xa space is not obtained simply by "neglecting" the other 
variables xi, in /prior, i-e. by setting Xb to the observed values in /prior: 

)■ (47) 

Since the observer region O can in general have a complicated shape in Xi space, the effective 
distribution function can have a quite different form than the expression in the right-hand-side 
of Eq. (47) . The effective probabihty force defined in Xa space is then given by 

Feff = V/eff, (48) 

where V = {d/dxi, • • • , d/dxn). 

The argument here suggests that it is a rather common phenomenon to have a nontrivial 
probability force when we focus only on a subset of the entire parameter set Xj, which is almost 
always the case. Then, if the resulting probability force is strong, we are likely to encounter an 
observer naturalness problem in the form of (b) or (c) in Fig. 7. 

6.3 Cut factor and comparisons between different theories 

In general there are many theories T that lead to the low energy Lagrangian with parameters x 
and that are described by a multiverse distribution frix). As an example, x may be the Yukawa 

Depending on the shape of O, we may have to integrate xn over multiple domains for some values of 
xi,--- ,xn-i- 
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couplings of the up and down quarks and electron, yu,d,e-, a-nd T may label the various theories of 
flavor. The number of observers in the multiverse who are governed by theory T is 

riT = j fxix) n{x) dx. (49) 

A typical observer will be governed by the theory T that has the maximum value of ut (assuming 
that several T do not have comparable ut)- The absolute normalization of n{x) is unimportant, 
since we are only interested in relative numbers of observers governed by different theories 

riT ^ lofTjx) n{x) dx ^ Nt_C^ 

^T> ~ la h' {x) n{x) dx ~ Mt' C^} ' ^ ' 

where Mt = J fri^) dx is the total "number" of universes (including volume factor weights) 
described by theory T and 

~ f fyfa;) dx , , 

Ct = — ~ . (51) 

JofT{x)n{x) dx 

Suppose now that n{x) is relatively fiat over O so that it is a good approximation to take it 
constant. Multiplying this constant by Ct then gives 

Cr = ikM^. (52) 
loMx) dx 

In this case there are a factor Ct more universes described by T in the multiverse than in O. 
Environmental selection solves observer naturalness problems at a cost of removing a factor Ct 
of these universes — we call Ct the cut factor for theory T. 

Consider, for example, the following two theories — the Standard Model where only the 
Higgs mass-squared parameter scans with a distribution function constant in |m^| and a theory 
beyond the Standard Model where only the weak scale scans with a distribution function constant 
in l/\nv. Considering selection by Eq. (29), the cut factors in these theories are 0(10^^) and 
O(IOO), respectively. If we take this naively, it seems to suggest that the latter theory is preferred 
over the former, and we might say that the Standard Model is less likely since it involves the 
cost of a much larger cut factor. (This result is similar to that of the conventional naturalness 
argument, although here observable universes always have v <^ Mpi due to selection effects.) Of 
course this is not rigorous — our ignorance of Mt, as well as other possible selection effects, 
makes the argument based only on the cut factor unreliable; without knowledge of A/r there is 



^^Here we do not consider the cut factor arising from selection of the cosmological constant, A <> Q T^^- In 
fact, the size of the cut factor associated with the cosmological constant does not depend, e.g., on the scale of 
supcrsymmctry breaking if the constant term in the supcrpotcntial, Wo, scans up to the fimdamcntal scale, M*, 
with the distribution flat in |Wop, as suggest by string theory [24]. This part of the cut factor is then of order 
Mf/A pa 10^^° regardless of the theory. (An exception for this is given by supersymmetric theories in which both 
supersymmetry and R symmetry are broken at low energies.) 
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no real reason to prefer a theory beyond the Standard Model over the Standard Model (other 
than considerations based on other physics such as gauge coupling unification and dark matter). 
Nevertheless, the relative cut factors for two theories do contribute to the relative number of 
observers described by each theory, as in Eq. (50), hence we may describe a large cut factor for 
some theory as a cost in its solution to an observer naturalness problem. 

6.4 Evidence for environmental selection 

We have seen that observer naturalness problems, which appear in the form of (a), (b) or (c) of 
Fig. 7, can be solved in general by environmental selection. Turning this argument around, we 
find that effects of environmental selection can show up in one or both of the following forms: 

(i) Amazing coincidences: we are living in a region of parameter space that admits complex 
observers, which, however, is only a (very) small portion of the entire parameter space. 

(ii) Living on the edge: we are living (very) near an "edge" of the parameter region O that 
admits complex observers, i.e. physical parameters take values that correspond to a point 
very close to the observer boundary, compared with the extent of the region O. 

In fact, how environmental selection manifests itself depends on the basis — we have seen that 
a phenomenon that appears as (ii) in one basis can appear as (i) in another basis, or vice versa. 
However, since we are often presented with a natural basis in which the physical meaning of 
parameters is most direct and/or intuitive, it is meaningful to consider the manifestation of envi- 
ronmental selection in that particular basis (physical basis). As we have seen, the manifestation 
then takes the form of (i) if the effective distribution function fesix) is (nearly) constant over the 
observer region O, while it takes the form of (ii) or a combination of (i) and (ii) if fesix) is peaked 
towards a boundary of O. There are many possible origins for a nontrivial form of fef[{x): an 
X dependence of the landscape vacua, the dynamics of the population mechanism, the observer 
factor n(x), and environmental selection acting on variables other than x. 

Can the observation of one or both of the above two phenomena, (i) and (ii), be viewed as 
evidence for environmental selection? In a given theory, T, the answer is yes, since otherwise it 
is very difficult to explain these features. Even though there is an observer naturalness problem 
associated with a very small value for Pt, after environmental selection a typical point within O 
results, i.e. 

Pt Po,T ~ 1- (53) 

The smaller the original Pt-, the more amazing the coincidence, or the closer to the edge we are 
living. The significance for the evidence for environmental selection is then quantified by the size 
of the naturalness probability Py, with a smaller value of Pt corresponding to stronger evidence. 
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In practice, however, we do not know beforehand the correct theory describing our universe. 
Suppose we encounter a situation described by (i) or (ii) above. Do we conclude that we have 
found evidence for environmental selection, or that the theory we are considering is simply wrong? 
Whether an observation of the phenomenon described in (i) or (ii) - an observer naturalness 
problem - can be viewed as evidence for environmental selection depends on whether one can 
find an alternative theory in which the problem docs not arise. A simple alternative theory without 
the naturalness problem may provide a better description of our universe. On the other hand, it 
is possible that we cannot find such a theory, or can find only theories that are significantly more 
complicated. Then we may conclude that environmental selection provides the best explanation 
of the phenomenon, so that we provisionally accept the multiverse theory. 

In principle, with competing theories to describe nature, the evidence for environmental se- 
lection at a particular observer boundary is given by the largest value of Pt associated with that 
boundary. However, it may be that a theory with large Pt is very complicated or ad hoc, so that 
it does not appear to be an adequate solution to the observer naturalness problem. In this case 
it may be judged that the evidence for environmental selection is better represented by a smaller 
value of Pt associated with some simpler theory. 

The discussion given above illustrates why the cosmological constant is a powerful argument 
for environmental selection: in all known quantum field theories the naturalness probability is 
extremely small (Pt ~ 10~^^° - lO"^'^ depending on the presence of weak scale supersymmetry 
and the nature of its breaking), implying that the largest value of Pt is of order 10~^° or smaller. 
The absence of a simple theory with Pt — 0{1) is crucial in this argument. A similar argument 
may also be made for the quark and lepton masses discussed in sections 4 and 5. We have argued 
that all known theories of flavor arc quite inadequate to explain the relevant observer naturalness 
problems, leading to Pt ^ (10^'^ - 10^^). Although this is numerically less impressive than the 
case of the cosmological constant, it is nevertheless very important. A single piece of evidence for 
environmental selection, no matter how signiflcant, could be completely erased by the discovery 
of a single new theory The more independent pieces of evidence for environmental selection, the 
more convincing the overall picture becomes. 

One might think it always difficult to "confirm" the absence of alternative theories which do 
not have the problem. Indeed, the absence of such theories can in general be inferred only from 
negative results of theoretical search. However, in the case that the observer naturalness problem 
is related to a fine-tuning naturalness problem, it is possible that we can be convinced rather 
firmly that the problem does in fact exist. This will be the case, for example, if we do not see 
any deviation of gravity from Newton's law down to a scale (much) smaller than 0(100 //m), 
since it will tell us the absence of a physical threshold that can control the observed value of 
the cosmological constant. In sections 8 and 9, we will also argue that the observation (or non- 
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observation) of new physics at the TeV scale may also be viewed, depending on what we will 
see, as evidence for the existence of an observer naturalness problem, and hence environmental 
selection. 

We stress that environmental selection can provide numerical predictions that arc difficult to 
obtain in other ways. Specifically, this occurs if the effective distribution function fcs has a non- 
trivial form in the physical basis. In this case the physical parameters take values corresponding 
to a point close to the observer boundary, giving nontrivial predictions. 

An example of predictions made possible by environmental selection is given by the stability 
boundary of the desired electroweak phase of the Standard Model [15]. Suppose that the Standard 
Model is valid up to some high scale near the Planck scale and that the weak scale v results 
from environmental selection. Suppose further that the Higgs quartic coupling at M*, A/^,*, 
and the top quark Yukawa coupling at M*, vary from one universe to another. There is 
then an observer boundary 0{\h,*,yt,*) — 1 corresponding to sufficient stability of the desired 
electroweak phase. Now, if the distribution function fesi^h,*,yt,*) is strongly peaked towards 
the phase boundary then our universe is expected to be close to this edge. In particular, if the 
peaking is stronger in A^^* than in * then the most probable point on the phase boundary has 
-^Higgs — 107 GeV and nit — 175 GeV. Discovery of the Higgs boson near this mass, together with 
the absence of any new electroweak physics, would then provide evidence that our universe is near 
the edge of this observer phase boundary, and hence of environmental selection. In general, if we 
find ourselves living close to an observer boundary and if we do not have a simple (alternative) 
theory explaining that fact, then we may regard it as evidence for environmental selection. 

7 Predictions for m„, and rrie from Environmental 
Selection 

The observer naturalness problem associated with the stability of neutrons, deuterons and com- 
plex nuclei was introduced in section 4, and is illustrated in Fig. 5. We showed that, no matter 
what the theory of flavor, there is always a factor in the naturalness probability, P, from the depen- 
dence on f/AqcD- In section 5 we argued that there is a factor, Pp, in the naturalness probability 
that is highly dependent on the theory of flavor. In all known theories P = PpP ^ (10^'^ - 10~^). 
In this section we argue that this naturalness problem can be solved by environmental selection; 
furthermore there arc several possible solutions with different consequences. 

We begin by restating the naturalness problem in terms of Standard Model parameters. For 
the neutron and complex nuclei stability boundaries a crucial quantity is 

m„ - mp - me = Ci{yd - yu)v - yeV - CaaAgcD, (54) 
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Figure 11: Q value for the reaction n ^ pe u. 



where Cj ~ (0.5 - 2) and Ca ~ (0.5 - 2) are strong interaction coefficients, and we choose 
a definition of the QCD scale such that Aqcd = 100 MeV. Throughout we use approximate 
ranges for m„ o and md^o from Ref. [20]. For complex nuclei the binding energy per nucleon is 
-E'bin = CfiAqcD, with Cb ~ 0.08 another strong interaction coefficient. The neutron and complex 
nuclei boundaries are then described by the inequalities 

V V 

< Ci{yd - Vu)-. 2/e-r C'^a < Cb- (55) 

^^QCD Aqcd 

In the approximation that the deuteron binding energy, Bo-, is linear in m„ + nid in the region of 
interest, the stability boundary for the deuteron can be written as 

= Ci - C^iVd + Vu)^ > 0, (56) 
Aqcd Aqcd 

where Ci^2 are two further strong interaction coefficients. These three boundaries involve just 
four independent combinations of Standard Model parameters, yu,d,eV / Aqcd and a. 

Much of the observer naturalness problem arises because the three terms that contribute to 
TTLn — rrip — rrie in Eq. (54) are comparable, as shown in Fig. 11, even though Yukawa couplings 
and the ratio of mass scales v/Aqcd are expected to range over many orders of magnitude. 
The electron mass, i/eV ^ 0.5 MeV and the electromagnetic mass splitting 6em = C'a«AQCD — 
1.0 =F 0.5 MeV differ by only about a factor of 2. They both stabilize the neutron, and hence 
an unstable neutron requires i/d > Vu- Not only is this the case, but the isospin breaking term 
5d-u = Ciirud — rriu) = 2.3 ± 0.5 MeV only just over-compensates the negative terms to give a 
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net available energy for neutron beta decay of just 0.8 MeV, an amount that is also very close to 
the electron mass and electromagnetic terms. In any theory of flavor where yu,d,e are determined 
by symmetries, these numerical coincidences, allowing neutron instability and the existence of 
complex nuclei, must simply be accidents. Furthermore, these accidents also involves v/Aqcd and 
a. The other Standard Model combination, {y^ + Vujv/ Aqcb, which could also range over many 
orders of magnitude, is numerically within about a factor of 2 of the deuteron stability boundary. 
Finally, as is apparent from Fig. 5, the observed parameters arc actually much closer to the 
neutron stability boundary than a factor of 2. The perpendicular distance in coupling parameter 
space is (10 - 30)%, depending on the values of Cj and yu,o/yd,o- As discussed in previous sections, 
some aspects of these accidents could arise from a theory of flavor that determines yu,d,e to be 
small and in the approximate ratio 1 : 2 : 1 at a unified scale; but the accidents associated with 
the value of v/Aqcd and the closeness to the neutron stability boundary remain. 

Environmental selection cuts out large region of parameter space where neutrons are stable, 
or deuterons or complex nuclei are unstable, greatly reducing the observer naturalness problem. 
The only remaining question is whether our observed universe is a typical universe within the 
observer region O shown in Fig. 5. This depends on the probability distribution for yu,d,eV / -^qcb 
and a, and on whether O has been correctly identified, allowing several solutions of the problem 
with different implications. We begin by taking a fixed to its observed value while allowing yu,d,e, 
V and Aqcd to scan. Since there are only two scales v and Aqcd appearing in the nuclear stability 
observer boundaries, we can take Aqcd = 100 McV without loss of generality, setting the unit of 
mass. Furthermore, since v always multiplies a Yukawa coupling, environmental selection for the 
light fermion masses can be discussed in terms of a distribution function f{mu,d,e) without loss 
of generality. ""^^ For a given / we can compute (m„), (m^i) and (mg) and compare them with the 
values observed in our universe, m„ o, m^j „ and rrie^o- We can also compute how close a typical 
universe in O is to the neutron stability boundary. 

Predictions for {niu^^e) involve the mass scales that arise in the observer boundaries, namely 
the maximum value of m+ = rud + m„ allowed by deuteron stability 

m+^ax = ^Aqcd ^ (1.4 - 2.7) m+,„ ~ (1.4 - 2.7) x (5 - 12) MeV, (57) 

for a — (1.3 - 5.5) MeV, the binding energy per nucleon in stable complex nuclei 

^bin = CbAqcd ~ 8 MeV, (58) 

^^The distribution function f{mu,d,e) here and below really means f{mu,d,e/-^QGu), but we omit Aqcd for 
notational simplicity. In the language of section 6.2, this is the effective distribution for mu,d,e/^QCD after 
integrating out the other parameters. Specifically, /(m„,d,e/AQCD) = (5(m„/AQCD - tJuv/Aqcd) 5{md/ Aqcd - 
Udv/ Aqcd) S{me/ Aqcd - yev/AqcD) f{yu, yd, Ve, v, Aqcd) dyu dyd dye dv cJAqcd- Using /(m„,d,e), we are able to 
discuss environmental selection at nuclear boundaries without answering how the degeneracy inside mu,d,e (scaling 
u/Aqcd and yu,d,e oppositely keeping mu,d,e fixed) is determined. 
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and the electromagnetic contribution to the proton mass 

Sem = CaaAQCB ^ (0.5 - 1.5) MeV. (59) 

While all three mass scales are proportional to Aqcd, ^em is significantly smaller than m+max and 
E'bin- Hence it is important to see which of these mass scales enter the predictions for {rriu^d^e). 
We consider three situations that solve the observer naturalness problem: 

I. An important part of the observer boundary is missing. 

II. The probability distribution is flat in mass space, df/ = Adm^^drnddme- In this case the 

closeness to the neutron stability boundary discussed above is accidental. 

III. The probability distribution dAf = f{mu,md,me) dmudmddme yields a nontrivial probabil- 
ity force towards the neutron stability boundary. 

For case I, consider a multiverse with m^, and irie uniformly distributed on logarithmic 
scales so that the relevant diagrams are shown in the left panels of Fig. 5. Although the original 
naturalness problem P ^ 1 has been ameliorated by a large cut factor, the naturalness problem is 
not entirely removed by environmental selection for neutron instability and deuteron and complex 
nuclei stability, since Pq ^ 1. For example, there arc large regions of O at small m„ e that are 
distant from the observer boundary. A complete solution follows if there arc additional relevant 
boundaries which we have failed to identify, that reduce O to the point where our universe 
becomes typical within O. This would certainly require new physical constraints to remove the 
large regions with low values of niu/mu^o and me/me,o- A complete solution may need further cuts 
to remove large values of me/rrie^o not already excluded, for example using the threshold for the 
pp reaction. Furthermore, to understand our closeness to the neutron stability boundary, it would 
be necessary for other cuts to approach our universe with a similar closeness on a logarithmic 
scale. While we can certainly identify physical processes that would introduce extra boundaries, 
we are unable to argue that they induce catastrophic changes rather than just the substantial 
changes discussed in section 4. 

For the second case, II, suppose that the distributions of m^, and nie are flat on linear 
scales, so that the effects of environmental selection can be understood from the cuts of the 
observer boundaries drawn in the right panels of Fig. 5. The observed masses are relatively 
typical within O, so that Pq is not much smaller than unity, and the naturalness problem is 
largely solved. In fact, the observed masses are still close to the neutron stability boundary, even 
in the right panels of Fig. 5, which in this example is accidental. Having assumed a simple form 
for the multiverse distribution, i.e. that /(m„,md,me) is constant, we are able to use the precise 
form of the observer boundary to compute the average observed values of the electron, up quark 
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and down quark masses by integrating over O 

(me) = \C: m+^a. ~ C, (^^) (2 - 3) MeV, (60) 
3 

(m+) = -m+max ~ (1 - 2) m+,o, (61) 
1 

(m_) = -m+max ~ (0.7 - 1.4) m+,o, (62) 

where m± = md±mu. In these equations the analytic expressions are obtained without including 
the effect of the complex nuclei boundary, and the numerical range corresponds to a = (1.3 - 
5.5) MeV. (Including the effect of the complex nuclei boundary changes the numerical values only 
up to about 30%.) These results demonstrate that environmental selection yields predictions for 
mu,d,e once a simple form for the probability distribution has been assumed. The predictions for 
rUu and are good. The prediction for nie is quite uncertain. For many values of the strong 
interaction parameters it is somewhat large; for example, for central values of C/, a and m+,o, 
(me) ~ 6me,o, so that low values of C/ and m+,o as well as a high value of a are preferred. Never- 
theless, we stress that a major part of the observer naturalness problem is solved if f{mu, m^, me) 
is relatively fiat within O. On linear scales for mu,d,e, our universe is quite typical of O. Of course, 
case II does imply that the closeness to the neutron boundary is accidental, and the rest of this 
section is devoted to understanding this closeness. A peaked distribution function can also give 
(me) oc 5-EM, rather than oc m+max, leading immediately to an understanding of the lightness of 
the electron. 

For case III, inside O the distribution / is peaked towards the neutron stability boundary, 
allowing us to explore the consequences of an environmental explanation for why mu,d,e are so 
close to this boundary. We call this stability boundary the n surface — it is a 2-dimensional 
plane in the 3-dimensional space of masses, mu,d,e- Within O, the probability force F = V/ can 
be resolved into components parallel, F\\, and perpendicular, Fj_, to the n surface. We assume 
that Fj_ points towards the n surface. The Fy field will determine the most probable location on 
the n surface. 

The n surface has a triangular shape, as shown in Fig. 12. Two of the edges of the triangle 
correspond to edges of physical space, mu,e = 0, while the other edge corresponds to the inter- 
section of the n surface with the D surface, the boundary for deuteron stability. We call these 
three edges the u, e and D edges, and label the three vertices of the triangle as ue, uD and eD, 
as shown in Fig. 12. The Fy field will determine where on the n surface triangle the distribution 
/ is maximized; there are three possibilities: 

III-l. In the interior, not close to an edge. 

III-2. On an edge, not close to a vertex. 
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Figure 12: 2-dimensional n surface in the 3- dimensional space of masses, mu,d,e- 



III-3. On a vertex. 

On a linear scale, the n surface is small compared to the expected available parameter space, so 
that it is unlikely that / would have a sharp peak in the interior of the triangle. However, the 
observer factor, n, may vary significantly across the triangle inducing a maximum in the interior. 
For case III-2, as we move along an edge, / rises to reach a maximum and then falls. As an 
example consider the distribution in O to be 

f{mu,md,m^)=A{md-mu)~^, (63) 

where yl is a normalization constant and p is positive. This situation would arise if the populated 
landscape has a low probability to yield a universe with large breaking of isospin, or if there 
are more observers in a universe as isospin is restored. In a slice through parameter space at 
constant rrie, the resulting probabihty force, F_ = V/, is towards and perpendicular to the n 
surface. However, in the full 3-dimensional space F_ is not perpendicular to the n surface. This 
is apparent in Fig. 13 which shows the n surface in the 3-dimensional space of m+, m_ and 
nie- Clearly the force will lead to a preference for low values of rrie, so that if p is large enough 
the most probable universes will have mg close to zero, i.e. close to the e edge of the n surface 
triangle of Fig. 12. From Fig. 13, all points on the e edge are equally probable, so that there is 
no expectation of being close to either the ue or eD vertices. We can compute the expectation 
values of mu,d,e by integrating over the 3-dimensional region O. For p > 3 the regions in Fig. 13 
at low m_ and low rrie dominate the integrals, and one sees that 6em sets the scale for both rrie 
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Figure 13: The n surface in the 3- dimensional space of m+, m_ and mg (shaded by sohd hues). 
The iTLe axis is stretched relative to the m± axes to make the figure more visible. The observer 
region O is the region surrounded by the n surface, the m„ = surface (shaded by dashed lines), 
the nif. = surface, and the deuteron stability surface at = m_|_jnax- The complex nuclear 
stability boundary is not shown. 



and m_, giving 

(me) = Sem, (64) 

p — 3 

and ^ 

Ci (m_) = - — - 6em- (65) 
p — 6 

This explains why the three contributions to m„ — rrip — rrie shown in Eq. (54) are comparable. 
On the other hand. Fig. 13 shows that m+ is uniformly distributed along the n surface from small 
values to m+max, and hence 

(m+) = - m+max- (66) 

Notice that, of the three mass scales that enter the observer boundaries, Eqs. (57, 58, 59), i^^bin 
does not appear. This is because the probability distribution is peaked away from the complex 
nuclei boundary, which therefore becomes irrelevant in determining the averages. 

For case III-3 there are three possible vertices of the n surface triangle to consider. Our 
universe is certainly not close to the uD vertex, since this gives rrie too large, hence we study the 
remaining two vertices. There are many / that lead to these vertices; we begin by identifying 
and studying the simplest cases. The ue vertex can be reached by a preference for low values of 



^"^The appearance of the mass scales in Eqs. (64, 65, 66) depends only on a probability force towards low m_, 
and not on the particular choice of the power law / in Eq. (63). 
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n^u,d,e- However, low values of m„ or rrie are not sufficient as they lead to the u or e edges, hence 
as a simple example we study the distribution 

/(m„, md, me) = A ^ (67) 

with q positive. At the ue vertex the value of is determined by 5em- For sufficiently large q 

this is the only relevant scale so that the average values of niu^d^e over O will all be determined 

by Sem- ^ 

(me) = Sem, (68) 

g — 4 

Ci (m+) = — ^ 5em, (69) 

C,(m_) = ^(5em. (70) 

While choices of q and Cj can lead to a hierarchy between these averages (for example (me) ~ Sem, 
(m+) ~ 10 5em and (m_) ~ 6 5em for g = 5 and C/ = 0.5), the hierarchy is more readily generated 
when (me) oc 5em and (m+) oc m+max, as in Eqs. (64, 66). 

The eD vertex is favored by large m^ and small me] however large or small m^ both lead away 
from the eD vertex. A simple distribution peaking towards the eD vertex is 

/(m„, ma, me) = A m[^, (71) 

with r ^ m+max/f^EM- Any distribution peaking near this vertex will lead to (m_) oc Sem and 
(m+) ~ m+max, which agree with observations. The prediction for the electron mass depends on 
the strength and direction of the force. In the example of Eq. (71), (me) ~ m+max/'"- 

In two Higgs doublet theories a probability distribution for the ratio of vacuum expectation 
values, tan/3, contributes to the distributions for mu,d,e- Suppose that this distribution favors 
large tan/3 and that tan/3o 3> 1. In this case low nid and mg are preferred, so that within the 
observer region the probability distribution is peaked towards the ue vertex of the n surface, 
leading to (me), (m+), (m_) oc 5em- The probability distribution in the electroweak symmetry 
breaking sector may be the origin of the force determining {rriu^d.e)- 

Finally we consider the possibility that a also scans. In this case the size of the electromagnetic 
mass difference, 6em, scans relative to the purely QCD scales of E'bin and m+max- This means 
that there is a shift in the position of the allowed window for (C/m_ — me) from neutron and 
complex nuclei stability. If a increases too much this window shifts to a non-physical region where 
m_ > m+, leading to an upper limit on a 

a < (72) 
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which is numerically about 0.2, with large uncertainties. Thus a is about an order of magnitude 
away from the maximum value that it may take anywhere in O. 

If a is very small then the range for [Cirri- — rrig) increases. The amount of increase is 
negligible on a linear scale, but is sizable on a logarithmic scale. For example, with a = 10~^, 
{Cirri- — rrif.) can range from 10~^ MeV to 8 MeV. However, if gauge couplings unify then a and 
Aqcd become related. As a is decreased, so Aqcd becomes exponentially smaller. It could be 
that selection effects on the size of Aqcd compared with the unified scale and/or the electroweak 
scale dominate over selection effects of a in nuclear physics. 

In many circumstances, precise predictions for environmental selection follow from assuming 
sharply varying distribution functions /. However, in the present example of nuclear physics 
in the parameter space mu,d,e and a, it is possible that / is sufficiently slowly varying over O 
that physical arguments based on the observer factor, n, could lead to predictions without any 
assumptions of sharply varying /. For example, consider variations in the parameters within 
O that lead from the neutron surface to the observer boundary with no complex stable nuclei. 
Moving towards the boundary of no complex stable nuclei, the observer factor may be reduced 
by successive nuclei becoming unstable. On the other hand, moving closer to the neutron surface 
than our universe will lead to a longer neutron lifetime and therefore to more primordial helium 
production; the reduction in primordial hydrogen will result in fewer hydrogen burning stars. The 
competition of these (and other) effects may determine the location of our universe within O. Of 
course, even in this case, some assumption about / is still needed. 

8 Electroweak Symmetry Breaking selected by Nuclear 
Stability 

The origin of electroweak symmetry breaking is one of the largest mysteries remaining in the 
Standard Model. The quadratic divergence of the Higgs mass-squared parameter in the Standard 
Model implies that if the scale of new physics M is (much) larger than the weak scale v, the theory 
requires fine-tuning. This hierarchy problem was a key motivation for much of the model building 
in the last 30 years. What do we know about the scale M? In many (non-supersymmetric) theories 
beyond the Standard Model, precision electroweak data indicates a "little hierarchy problem": 
v/M is uncomfortably small, typically of order (10~^ - 10~^) or smaller (sec e.g. [25]). There 
is also a similar fine-tuning problem in supersymmetric theories, although its origin is different. 
A sufficiently heavy Higgs boson typically requires superparticles to be somewhat heavier than 
the weak scale, leading to some amount of fine-tuning (see e.g. [26]). In either case, we find 
that some amount of unnaturalness is present, at least for the simplest theories, suggesting that 
environmental selection may be playing a role. In this section we investigate whether a hierarchy 
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between v and M is to be expected from environmental selection, the size of any such hierarchy, 
and how the hierarchy depends on which parameters are assumed to scan. 

To address the question of the environmental selection of v/M, several issues must be ad- 
dressed: what is the theory under discussion, which parameters of that theory scan, and what 
observer boundaries implement the selection? Below we formulate a fairly general class of theories 
that describes clcctrowcak symmetry breaking, and the nuclear stability boundaries of section 4 
arc used to implement selection. In section 9 we consider an alternative possibility that selection 
occurs at the clcctroweak phase boundary. 

The mass scale of the new physics that generates electroweak symmetry breaking is defined 
to be M, and we assume that the effective theory below M is the Standard Model with the Higgs 
potential 

V = mlh^h+^{h^hf. (73) 

Integrating out the physics of the electroweak symmetry breaking sector at scale M in general 
leads to several contributions to ml, some positive and some negative, which we write as 

ml^{g^{xi)-g2{xi))M^, (74) 

where the functions gi_2 are both positive. The dimensionless parameters ,Xj are the set of pa- 
rameters of the theory above M that substantially affect electroweak symmetry breaking. These 
parameters are evaluated at the scale M, so that gi_2 arc also functions of M through the logarith- 
mic sensitivity of Xi on M: gi 2{xi{M)). Note that here includes the quadratically divergent 
radiative corrections in the Standard Model that are regulated by the theory above M. The 
parameters Xi thus include the Standard Model gauge and Yukawa couplings. 

The precise nature of the couplings Xi and the functional form of gi^2 are unimportant for our 
discussion, so that we write gi{xi) = Ax and g2{xi) = Ay, giving 

ml^{x-y)AM', (75) 

where x,y,A > 0. The numerical constant A is chosen such that typical values for y in the 
observer region are of order unity. The parameters x and y depend logarithmically on M through 
renormalization group evolution. {A is a one-loop factor in many theories beyond the Standard 
Model. Since ml contains quadratically divergent contributions in the Standard Model, and 
thus X and y contain pieces proportional to the SU{2) gauge and top Yukawa couplings squared, 
respectively, the value of A should not be much smaller than the one-loop factor.) In the context 
of conventional naturalness criteria, electroweak symmetry breaking is unnatural if x is typically 
larger than y in an ensemble (only a small fraction of members in the ensemble leads to electroweak 
symmetry breaking) , whereas if x is typically of order y or less then the natural value of the weak 
scale is y/A M for a quartic coupling of order unity. 
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Figure 14: The (supersymmetric) little hierarchy problem as an observer naturalness problem. 
Contours of M/v are drawn in a 2-dimensional slice of parameter space of a generic electroweak 
symmetry breaking sector. As the experimental limit on M increases, so the allowed region of 
parameter space shrinks to values of larger M/v close to the phase boundary. 



A key observation is that as experimental limits on physics beyond the Standard Model get 
stronger so the mass scale M is constrained to be larger. In some (non-supersymmetric) theories 
of electroweak symmetry breaking, this typically arises from contributions of particles of mass 
M to the precision electroweak observables. In many supersymmetric models, the increased 
lower limits on the Higgs boson and superparticle masses have pushed up the mass scale of some 
superparticles to about a TeV or larger. These little hierarchies between v and M imply that the 
parameters Xi are constrained towards the phase boundary of electroweak symmetry breaking. 
As experiments push up M, so gi and g2 of Eq. (74) cancel to give v/ \/~AM (much) smaller than 
unity. This is illustrated in Fig. 14 for a 2-dimensional slice through the parameter space. The 
unusual closeness of the parameters to this boundary can be viewed as an observer naturalness 
problem. 

In this section, we investigate the various possible ways in which environmental selection on 
a multiverse can play a role in electroweak symmetry breaking. We use the class of theories 
described above, which contains most non-supersymmetric and supersymmetric theories beyond 
the Standard Model. The relevant observer boundaries are those of nuclear stability, shown in 
Fig. 15, so that in general one must consider scanning parameters that affect nuclear physics as 
well as the parameters of the electroweak symmetry breaking sector. We start by considering 
only a few parameters scanning, and then progress to more general situations. 

First of all, it is natural to expect that some parameters of the theory at M scan in the 
multiverse, so that there are universes with different values of the weak scale v. Universes with 



49 



0.5 1 1.5 2 2.5 

mu/nriu.o 



Figure 15: The observer region in m„-md space (inside the shading). The dashed hne drawn from 
the lower left to the upper right represents the trajectory followed when only the electroweak 
vacuum expectation value v is varied (neglecting the effects of a variation of me) . 

X > y have v = 0, while those with y > x have 

Scanning of parameters in the electroweak symmetry breaking sector at M may also lead to 
a nontrivial effective distribution for Xh{v). Since a light Higgs boson has not been discovered, 
however, this distribution should lead to Xh{v) typically having an order unity value in the observer 
region. This implies 

in our universe, where y is of order unity due to definition of A. The crucial question is if 
some cancellation between y and x is expected, which cannot be understood in the conventional 
symmetry approach. If not, then v ~ y/AM as in the conventional case. If so, however, we obtain 
an extra hierarchy between v and M that cannot be explained by a symmetry, and the question 
becomes: is y — x typically of order 10~^-10~^, with the multiverse generating a little hierarchy, 
or is it extremely small, for example 10~^°, giving a large hierarchy? 

In addition to v, the nuclear stability observer boundaries depend on parameters yu,d,e, oc and 
Aqcd- In sections 8.1 and 8.2 we assume that these additional parameters do not scan, so that the 
only scanning in the Standard Model at low energies relevant for electroweak symmetry breaking 
is that of V. In this case, environmental selection for neutron instability and deuteron stability 
defines a fixed observer window for v 

<v < v+, (78) 
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where, from Eq. (29), f_ ~ 0.5 and u+ ~ 2vo. (The value of v+ depends on the parameter a 
describing the strength of the deuteron binding; here we take a — 2.2 MeV.) This is ihustrated 
in the mu-m^ plane in Fig 15. As v is varied about Vo so m„ and vary, but with a fixed ratio, 
as shown by the dashed line. (The corresponding variation of me gives only small effects.) Since 
Vo lies centrally in the observer window, it is consistent with a distribution /.„(';-;) that is slowly 
varying. Can environmental selection for v generate a little or large hierarchy, and if so is the 
effective distribution for v mildly varying? 

Another possibility, studied in section 8.3, is that the entire set yu,d,ei ol^ Aqcd and v scans, 
so that there is a non-zero probability distribution throughout the nuclear observer region. A 
2-dimensional slice through the observer region is shown in Fig 15; the scanning is no longer 
restricted to the dashed fine. An interesting question then is whether the closeness to the neutron 
stability boundary is a statistical accident, or whether it results with high probability due to 
a strongly varying distribution function. As stressed in section 4, the nuclear observer region 
depends on only four combinations of these quantities, mu,d,e/ ^q,CT> and a. For example, a 
common scanning of v and Aqcd does not affect nuclear physics. This implies that a numerical 
value for v is no longer determined by the nuclear observer region alone. Can environmental 
selection from the nuclear observer boundaries, shown in Fig 15, lead to a little or large hierarchy? 

8.1 Scanning the mass scale M 

The fundamental field theory at the cutoff scale M) will have a certain set of parameters. 

We assume that the scanning at is limited to those parameters that affect the scale M (which 
may be only M itself), leading to a distribution function f{M). In particular, the dimensionless 
parameters Xj(M*) do not scan. Assuming that effects on Xi from the parameters controlling 
M are small, this implies that the scanning of the parameters x and y in Eq. (75) comes only 
through a calculable logarithmic dependence on M, which arises from nontrivial scaling of these 
parameters under renormalization group evolution. 

Suppose now that y — x is positive and of order unity throughout the multiverse (which will 
be the case if y — x is positive at M^, and becomes larger as the renormalization scale is lowered). 
Environmental selection requires that v lies in the fixed range < v < v+. This selects M to be 
in the range 



There is, however, another possibility. It could be that as M varies so y — x passes through 
zero. In this case, there is a critical value M = Mc corresponding to the electroweak phase 
boundary, y — x = 0. By assumption, neither the critical value Mc, nor the nuclear window 




(79) 



The ratio between v and M is then the same as if there were no selection, v/M \/A. 
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Figure 16: Illustration of the ranges for M selected by the observer condition i>_ < v < v+. The 
solid curve gives v{M) in regions of M where electroweak symmetry is broken. 



parameters, v±, are scanning. It would be an accident for Mc to be close to v±, so we assume that 
they are distant. A physically interesting new possibility arises if — a; is small in our universe, 
which is possible if Mc ^ v±. This new possibility corresponds to the environmental selection of 
a large hierarchy. 

Such a large hierarchy can occur in two ways, depending on the sign of the beta function, 
f3, for illustrated in the two panels in Fig. 16. (The definition of (3 here is given by 

d{y — x)/d\nfi = (3.) For > (the left panel), electroweak symmetry breaking is only possible 
for universes with M > M^ and, since M^ ^ v±, the observer condition Eq. (78) selects M to be 
just above M^ in the range Mc(l + 5_) < M < Mc(l + 5+), where 6± = Xhvl/\p\AM^. The value 
of 1/ — a; is of order v'^/AM^ <^ 1, i.e. the cancellation of order v'^/AM^ is forced by environmental 
selection and the hierarchy between the weak scale and the scale of new physics is very large, 
Mc/v 1. (Note that x and y themselves are of order unity.) On the other hand, in the case 
that P < (the right panel) the broken phase has M < M^ and environmental selection gives 
Mc(l — 5+) < M < Mc(l — 5-); the value of y — x in this case is, again, of order v'^/AM^ ^ 1. 
For (3 > only a large hierarchy is possible. For /3 < 0, however, universes with v- < v < v+ 
are also possible with y — x ~ 1 and M ^ v/y/A, the case corresponding to conventional natural 
theories. The case of a large hierarchy is more probable if the distribution function for M, f{M) , 
is sufficiently weighted towards large M. 

One possibihty is that M is the overall scale of superparticle masses, so that the theory above 
M is supersymmetric. If the hierarchy is large, the situation described above then corresponds to 
the split supersymmetry scenario discussed in Refs. [7, 27]. However, the theory above M does 
not have to be supersymmetric. It may, for example, be a strongly interacting theory leading to 
a composite Higgs boson below M, possibly as a pseudo Nambu-Goldstone boson [28, 29, 30]. 
In this case we obtain a prediction on the Higgs boson mass as a function of M, by setting 
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Figure 17: Sketch of the 2-dimensional observer region in the M-z plane. 

the Higgs quartic couphng to be either very large or small at the scale M (compared with the 
logarithmically enhanced contribution between M and v). This could be interesting because for 
M <^ Mpi predicted values of the Higgs boson mass are outside the range of ^ (130 - 180) GcV, 
expected if the theory above v is the Standard Model up to a high scale of order the Planck scale 
with a (absolutely) stable electroweak symmetry breaking vacuum [31, 32]. 

We stress that the situation described here is very special: for example, if Aqcd scans in the 
multiverse, it is possible that the resulting hierarchy is "little," e.g. v/M pa (10~^ - 10~^), since 
the value of M may be set close to the observed weak scale by environmental selection, as we 
show below. 

8.2 Full scan of the electroweak symmetry breaking sector 

Next we consider a general scanning of the electroweak symmetry breaking sector, so that all 
of x,y,Xh and M vary independently, but we keep yu,d,e, ol and Aqcd fixed. In this case it is 
useful to consider a 2-dimensional observer region in the M-z plane, as shown in Fig. 17, where 
z = y — X. As before, we assume that \h{v) has a typical value of order unity, e.g. with its 
distribution function being strongly peaked at ~ 1, and hereafter we neglect the effect of its 
scanning. Electroweak symmetry breaking occurs in universes with z > 0, and the observer 
boundaries v — v± are shown in Fig. 17. Other parts of the observer boundary correspond to 

^^Our conclusions arc not affected by the scanning of Xh(v)- (The distribution of Xhiv) should, of course, be 
consistent with the bomid on the Higgs boson mass.) In some cases, for example in the minimal supersymmetric 
standard model, Aft(M) is determined by the theory at M and has very little ability to scan. 
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M — M*, the maximum value of M (the cutoff scale), and z = ^max, the maximum value of z 
determined by the ranges for x and y in the landscape, which we take to be of order unity. 

The effective probability distribution in this 2-dimensional observer region O is obtained from 
the distribution function for x, y, M as 

Mz,M)= I dxdyf{x,y,M)6{z-{y-x)). (80) 

With X and y scanning over the multiverse, the logarithmic dependence of these parameters on 
M frequently does not give a major effect, and hence we neglect it here. When these effects are 
important they lead to very interesting results, and we defer a discussion of this until section 8.5. 
If fcsiz, M) were constant over O, it is clear that our universe would be expected to have z ~ 1 
and V ~ ^/AM, since this corresponds to most of the area of O. However, it is also clear that, 
as fcs{z,M) becomes progressively more peaked towards low z and/or high M, so the typical 
hierarchy will first grow to a little hierarchy, z ^ 0(10^^ - 10^^), and finally to a large hierarchy, 
z -C O(10~^). A power distribution is easily able to overcome the narrowing of the observer 
region at large M. The probability distribution for the hierarchy, z, is given by 

Mz)= [ dMU{z,M), (81) 
Jo 

and a general scanning of the electroweak symmetry breaking sector of the theory could lead 
to a large probability for either a little hierarchy or a large hierarchy. Indeed, since a strong 
dependence of fcsiz, M) on z and M could result from many sources — the landscape of vacua, 
the population mechanism, integrating out other parameters and the observer distribution n — 
it would certainly not be surprising if either a little or large hierarchy resulted. Environmental 
selection of the weak scale can provide a sim,ple generic explanation for our present difficulties in 
constructing natural theories of electroweak symmetry breaking. 

Note that this mechanism works both in the context of non-supersymmetric and supcrsym- 
metric theories, and that the small value of z implies a cancellation between the positive and 
negative contributions in Eq. (74) {x and y) that cannot be explained in the conventional sym- 
metry viewpoint. In section 8.4, we study in detail the simplified case that y does not scan, and x 
and M have power- law or logarithmic distributions. We compute fz{z), find conditions for large 
and little hierarchies, and obtain an analytic result for the size of the hierarchy when it is little. 
We also compute the distribution for v^, fviv^), in the observer window. 

Before closing this subsection, let us consider the case that only the dimensionless variables 
in the electroweak symmetry breaking sector, with the value of M fixed. In this case 

experiments constrain the value of M to be larger than about a TeV. What value should we 

^^In supersymmetric theories there are at least two Higgs doublets, but our analysis can be applied to these 
cases by identifying h as the linear combination causing electroweak symmetry breaking. 
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expect for M? Without any special reason, we expect M to take some "random" value between 
TeV and fundamental scales; it is unlikely that M is close to the observed weak scale, since it 
requires an accident of order 0(0.01 - 0.1), as discussed in section 4.4. Hence wc typically expect 
a large hierarchy with M ^ v. The values of Xi are environmentally selected by Eq. (78) to a 
very small range: 

in which there is a large cancellation between x and y. 

8.3 Scanning over the entire nuclear stability observer region 

We continue to study theories of electrowcak symmetry breaking that lead to the Standard Model 
as an effective theory below M, with = {x — y)AM^. However, as well as having v scan, via 
the scanning of x, y and M, we now allow the other parameters of the nuclear stability observer 
boundaries to also scan. This means that instead of exploring a distribution through the observer 
region corresponding to only varying v, as shown by the dashed line of Fig. 15, we are now 
exploring a distribution over the entire observer region discussed in section 4. With yu,d,e, ol and 
Aqcd all scanning, the problem is apparently very complex. However, here we are only interested 
in two aspects of the problem: the probability distribution ior z — y — x that determines the 
size of the hierarchy between v and M, and the probability force perpendicular to the neutron 
surface that determines how close typical observers are to the neutron stability boundary. These 
questions can be addressed by studying an effective distribution over a reduced 2-dimensional 
projection of the parameter space. 

The neutron and complex nuclei stability boundaries of Eq. (55) can be written in the form 

i-<i< C+, (83) 

where 

A, ' 



f = (84) 



with 

M = £l^lhi_lA^M, (85) 
aAqcD 

and ^_ = Cq, and ^+ = Ca+CB/a- Note the similarity in the form of the equations for ^, Eqs. (84) 
and (83), to those for v, Eqs. (76) and (78). This means that the shape of the observer region 
in the 2-dimensional projection of the parameter space on the M-z plane, shown in Fig. 18, is 
the same as in the last subsection on the M-z plane. The curved observer boundaries at ^ = ^_ 
and ^ = however, now correspond to neutron and complex nuclei stability, respectively. Our 
universe thus lies very close to ^ — How probable this is will be determined by the effective 
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distribution for Similarly the expected size for the hierarchy will be determined by the 

distribution for z, fziz). 

In general, the probability distribution in our whole parameter space is given by the multi- 
dimensional distribution function f{x, y, M, X^, y^, y^, y^ Aqcd)) which is not easy to deal with. 
A crucial point, however, is that as far as questions of the size of the hierarchy, z, and the closeness 
to the neutron stability boundary are concerned, we only need to study the effective distribution 
function /efr(2;,M) obtained after integrating out all the other variables. In particular, we can 
parameterize our ignorance of the (potentially) complicated distribution function / and the shape 
of the observer region in this multi-dimensional space by the effective distribution function f^s in 
the 2-dimensional M-z space. Note that as long as the form of fes is kept arbitrary, there is no 
loss of generahty. 

The effective distribution function fes{z,M) is formally given by 

fes{z,M) = / dx dy dM dXh dyu dya dye da dAQCB f{x,y,M, Xh,yu,yd,ye, a, -^qcb) 
Jo 

xS{z-{y-x)) ^fM-^^^i^l^Mll^MV (86) 
V aAqcD / 

The problem of studying selection in the M-z plane then becomes identical to that in the M-z 
plane discussed in the previous subsection with the replacement M ^ M and v ^ ^, except that 
the upper boundary ^+ now depends on a scanning parameter a. The effect of the a scanning on 
the analysis, however, is small as long as the scanning is mild, e.g. the range of the a scanning 
does not span many orders of magnitudes. We assume this to be the case, and hereafter we 
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neglect the effect from this scanning. The effective distributions for z and ^ are given in terms 
of/eff(^,M)by 

/ dMMz,M), (87) 
Jo 

and 

f^{e)= jjzdM U{z,M)5(e - '-^), (88) 

respectively. 

We can identify three very different situations. The first is that fes{z,M) is very mildly 
varying. In this case, from Fig. 18 we see that a typical universe will lie in the middle of the 
observer region and hence environmental selection will lead to neither a hierarchy nor a closeness 
to the neutron boundary. A second situation has a strongly varying fcs{z,M), but with the M 
component of the probability force field unable to overcome the narrowing of the observer region 
at large M, shown in Fig. 18. In this case, if there is a strong force to low values of z, the 
multiverse yields a little hierarchy by making low values oi z ^ O(10~^ - 10"^) typical. If in 
addition there is a significant probability force to low M, the combined effects of the z and M 
distributions lead to a closeness to the neutron stability boundary as well as to a little hierarchy, 
as shown by the dots in Fig. 18. Finally, the probabihty force towards large M may be strong 
enough to give a large hierarchy with extremely small z. In this case the observer boundaries 
at ^_ and ^+ are very close to each other. At such low values of is it reasonable to have fes 
sufficiently different on these boundaries to favor the boundary at Here we must recall that 
/ is a product of a populated landscape distribution, /, and an observer distribution, n. It is 
certainly unreasonable for / to have such a large variation over such a small region of parameter 
space. However, as we move from the ^_ boundary to the ^+ boundary, nuclear physics changes 
very significantly, which is in a way independent of how small z is. Hence a closeness to the 
neutron boundary may be typical with a large hierarchy, but only if it is induced by n, through 
such arguments as appeared in the last paragraph of section 7. 

The closeness to the neutron boundary may be a hint that /cfj(-2, M) is not flat. As argued 
previously, there are many origins for strong probability forces, so that a little or large hierarchy 
is quite natural in the multiverse. In the next subsection we explore in some detail a subclass 
of the theories specified by Eq. (74), allowing explicit formulae for the effective distributions and 
the size of the hierarchy. 

8.4 An explicit example with power law distributions 

In the previous two subsections we assumed that the positive and negative terms in m^/M^ scan 
independently. To provide a simple explicit illustration of our results, in this subsection we choose 
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Figure 19: Sketch of the 2-dimensional observer region in the M-z plane, where z — 1 — x. 



to scan only the positive term, i.e. we fix y = 1 so that 

ml = {x- 1)AM'^, (89) 

and z — 1 — X. The analysis with x fixed and y scanning is very similar. For ease of calculation 
we assume a polynomial distribution function 

f{x,M)dxdM (xx^'dx MUlnM, (90) 

and the range of x is given by < x < Xmax, with Xmax some number larger than 1, while M has 
some very large maximum value M*. 

We begin by assuming that the only parameter appearing in the nuclear stability boundary 
that scans is t^, so that the condition of Eq. (78) gives the observer region 

< ^ {1 - x)^^ < vl- (91) 

The observer region O is sketched in the M-z plane in Fig. 19, which is obtained simply by setting 
y = 1 in Fig. 17. The region is separated by the nuclear physics observer boundaries from regions 
with V < V- and v > w^, and approaches very close to a; = 1, i.e. 2; = 0, as M approaches 

S> Vo- The narrowing of the observer region at large M is as we have seen in Fig. 17. 

What is the typical value of the scale Ml For a logarithmic distribution, g = 0, one might 
guess that all decades of M are equally probable; but this is not the case, because at large M 
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the observer region narrows. It is useful to compute an effective probability distribution for the 
hierarchy z — 1 — x — \mW/AM'^ over the observer region 

f,{z) = ^ dxdMfix, M)d{z-{1- x)) - (1 - zr (92) 

for sufficiently large M*. Prom the viewpoint of the hierarchy, this shows that the critical value 
of g is 2 not zero: when q — 2 the distribution for the hierarchy z is flat on logarithmic scales, 
at least for large M where z is small. Thus q > 2 gives a large hierarchy, shown by the narrow 
wedge at large M in Fig. 19, and q < 2 gives a httle (or no) hierarchy, corresponding to the bulk 
of the observer region in Fig. 19. For example, if the a priori distribution for M is logarithmic, as 
might happen if the scale is triggered by dimensional transmutation associated with some gauge 
dynamics, the most probable observed value of M is small because of the weighting from the 
narrowing of the observer region O}^ On the other hand, if dj\f ~ dM"^, as expected in a typical 
non-supersymmetric perturbative theory with mass scale M, we find g = 2, so that all scales are 
equally probable, implying a large hierarchy, but not one with M near M*. 

If g < 2, how large is the little hierarchy? For n > 0, the first factor in Eq. (92) suppresses the 
probability of having 2; ~ 1, i.e. the probability of having the lowest values of M gets suppressed, 
so that some amount of hierarchy between v and M arises from the scanning in the multiverse. 
For g < 2 the distribution of Eq. (92) leads to the average value of z: 



m 



2n + 4 - g 



^ (93) 



so that for 2n » 2 — g we obtain an extra hierarchy of a factor of ~ (2 — g)/2n from the 
multiverse, which cannot be explained in the conventional symmetry approach. Since |g| is 
expected to be small, this is an appreciable effect for large n. In Fig. 20 we plot fz{z) of Eq. (92) 
for (n, g) = (3,0), (10,0), (3, 1). We clearly see that the distribution is peaked towards z <^ 1. 
The probability of us observing fine-tuning of or more that cannot be understood by the 
symmetry approach is given by 

^ fn fz{z)dz , 

T^A = \i \' , (94) 
Jo Mz) dz 



which for A n ^ 1 reduces to 
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Va ^ ^, (95) 



^''In the case that M arises as a dimensional transmutation, the expected distribution of M below M* is not 
exactly cx d\nM, since it would lead to an unphysical conclusion that the probabilities of having M between 
10"M* and 10""'"^M, (n < —1) are equal for all n down to n ^ —00. (The distribution should, for example, be 
cut off at some value Mmin or in fact be a d{l/\nM); see discussions at the end of section 2.2 and in section 4.4.) 
The assumption/approximation of the exact logarithmic distribution, however, is sufScient for our purposes here. 
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Figure 20: The distribution function fz{z) for {n,q) = (3,0) (solid), (10,0) (dashed), and (3,1) 
(dotted). Each function is normahzed such that Jofziz) dz — 1. 




for q = (the case with dimensional transmutation) and 

for q = 1. For example, for n = 10 and g = (n = 10 and g = 1) we have a ~ 20% (~ 50%) 
probability of observing fine-tuning < 2% which cannot be explained in the conventional 
symmetry picture. 

The observer boundaries of Fig. 19 are the contours for v = so that, by construction, v can 
only be observed between f_ and What are the most probable values of v to be observed? 
The effective distribution for v"^ is defined by 

r (A A/f 2 \ 

U{v'') = J^dxdMf(x,M)5lv'-{l-x)^j. (97) 

For a large hierarchy with q >2, we find a flat distribution in 

fy{v'^) ~ const. (98) 

When yu,d,e, ci and Aqcd do not scan, the observed value of the weak scale Vo is roughly midway 
between V- and u^., as shown in Eq. (29), which is consistent with Eq. (98).^^ For a little (or no) 
hierarchy, with g < 2, we find a distribution 

Uiv") ~ v^-\ (99) 



-•^^Caution, however, is needed in interpreting Eq. (98). This result arises because at large M the contours of 
V- and va^ are extremely close in M-x space, so that the assumed form for /, Eq. (90), implies little variation 
between the contours. While this is expected for the multiverse distribution /, it may not be true for the observer 
distribution n. In this case Eq. (90) needs to be modified to incorporate the effect represented by n. 
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This is again consistent with Vg centrally located in its observer window, provided that q is not 
too negative. 

Next, we allow the other parameters that appear in the nuclear stability observer boundaries 
to also scan. As described in section 8.3, the issues of the size of the hierarchy and the closeness to 
the neutron surface can be addressed by studying a 2-dimensional projection of parameter space. 
Since y does not scan, z = 1 — x and so here wc discuss the M-x plane. We assume the effective 
probability distribution in this plane has a power law behavior 

fes{x,M)dxdM (xx^'dx M^dlnM. (100) 

The neutron and complex nuclei boundaries are again given by Eqs. (83, 84, 85), but with z — 1—x, 
and are represented in Fig. 18 in the M-x plane (with y now set to unity). The observer region 
as X — > 1 again becomes a narrow wedge where the hierarchy is large: M ':$> v. 

The equations describing this setup are very similar to those in the case with only v scanning, 
except that the observer boundaries are now at ^± rather than v±, and M is not the scale of the 
new physics, see Eqs. (83 - 85). The size of the hierarchy is governed by the effective distribution 
for z and is given by the analogue of Eq. (92): 

m-ii-zT^,. (101) 

Thus q > 2 gives a large hierarchy and q < 2 gives a little (or no) hierarchy. The size of the little 
hierarchy depends on n, and using the distribution Eq. (100) we find 

^ \ = (102) 

AM2 / 2n + A~q ^ ' 

The little hierarchy increases with n, as shown by the sequence of dots in Fig. 18 near the neutron 
stability boundary. 

What makes the proximity to the neutron boundary typical? The effective distribution for ^, 
which determines the probability force perpendicular to the neutron surface, follows immediately 
from a calculation analogous to that which led to Eqs. (98, 99) giving 

/^(^2-) ^ const. for q>2, (103) 

and 

~ e'-' for q < 2. (104) 

As with Eq. (98), caution is necessary in interpreting Eq. (103). At large M our assumed distribu- 
tion, Eq. (100), may not adequately account for variations in the observer factor n{^) between ^_ 
and In particular, we cannot conclude from Eq. (103) that a large hierarchy is incompatible 
with a multiverse explanation for the closeness to the neutron stability boundary. With a little 
hierarchy, sufficient closeness to the neutron stability boundary results with the distribution of 
Eq. (100) for -10 ^ g ^ -2. 
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Figure 21: Sketch of the 2-dimensional observer region in the M-z plane. The maximum value of 
-2) ^max(-^), is depicted by the dashed line. 

8.5 Landscapes with a critical value for M 

In the previous subsections we have ignored the logarithmic evolution of x and y. For many land- 
scapes this is permissible, but for some this evolution plays a critical role. Consider a landscape 
such that X* — y^, scans only in a restricted region with — > 0, where = x(M*) and 
= y{M^), which are the fundamental scanning parameters. In this case, electroweak symmetry 
breaking is only possible because of the evolution of x — y to lower energies according to the beta 
function j3 = d{x — y)/d\n.ix > 0. For each value of (x*,?/*), the quantity x(M) — y{M) passes 
through zero at some Mc(x*,|/*), so that in these universes electroweak symmetry breaking is 
only possible if M < Mc(x*,|/*)- As (x*,|/*) vary over the entire landscape, suppose that the 
largest value of Mc(x*,i/*) is Mc,maxj which is much less than M*. This maximum critical mass 
is clearly a property of the particular landscape under consideration; it defines a maximum pos- 
sible value for the electroweak scale anywhere in the multiverse, and has an important effect on 
the observer region of the electroweak symmetry breaking sector. For the critical universes with 
Mc{x^, y^) = Mc^max the evolution equation ior z — y — x can be solved as 

WM)=/3 1n^^, (105) 

where we have approximated that /? is constant. The trajectory of Eq. (105) is sketched in Fig. 21 
for the case that Mc^max ^ v±, which provides a new boundary to the observer region, since ^max 
represents the maximal value of z in the multiverse for a given M. In this subsection, we assume 
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that the distribution for M and z leads to most universes in the observer region having large M, 
close to Mc^max, SO that this new part of the observer boundary plays an important role.^^ 

If yu,d,ei and Aqcd do not scan then v± are fixed numbers that arise from nuclear selection. 
Since there is no reason for the critical mass Mc,max of the landscape to be close to a large 
hierarchy is expected. Here we are more interested in the case that v± varies in the landscape, 
and for simplicity we accomplish this via a scanning of Aqcd, which we take to have a mild 
distribution, for example flat on a logarithmic scale. In this case v± oc Aqcdi so that the observer 
boundaries in Fig. 21 labeled by v± have z± oc Aq^^j^, and hence move up and down as Aqcd 
scans. For a large hierarchy, M is so much larger than v that the interval of z that can contribute 
to universes in the observer region is very small Az ~ jM'^ . Since the distribution for Aqcd 
is mild, a gain in probability is accomplished by having the v± boundary curves of Fig. 21 move 
upwards, allowing a larger interval of z to contribute. As the v± curves move up so the values of 
M that are on the ^max observer boundary decrease. At some point this will imply a significant 
loss of probability via the distribution for M. For example, for f{M)dM oc M'^dlnM with large 
q, this loss of probability begins to set in once M ~ (1 — l/?)Mc,inax at which point z^^x ~ 
The mild distribution for Aqcd allows the v± curves of Fig. 21 to rise so that, at this value of 
M, z± will be close to z^ax.- A further drop in M causes a significant drop in probability from 
f{M), which is not offset by a sufficient growth in Az, since Zc{M) grows logarithmically. Hence 
we predict that a little hierarchy develops. Specifically, for a distribution function 

f{z,M,K(^CY,)dzdMdKqcYi^ f{z)M'idzd\TiMdhiKqcYi. (106) 

integrating over the observer region we obtain a typical size of the hierarchy 

{z) = ^'f^^^Jd' ^ ^, (107) 
ff{z)e-Tdz q 

where the range of the z integration is from to the maximum value that z takes in the energy 
interval between ^ v and M*, and we have taken f{z) to be constant in the last expression for 
illustration, although the result is not very sensitive to the form of f{z). Note that this result 
applies to an arbitrary positive value of g ^ 0{j3). 

Integrating over Aqcd effectively "removes" the v± boundaries in the 2-dimcnsional parameter 
space of {M, z}, and the resulting effective distribution in the M-z plane is unaltered for a prior 
distribution flat in In Aqcd- Allowing any mild prior distribution for Aqcd, or allowing other 
Standard Model parameters such as the Yukawa couplings yu,d,e to scan with mild distributions, 
would lead to an effective distribution fesiz, M) that would lead to the basic result of Eq. (107) 



^^For many landscapes the evolution oi x — y from M* to M can be ignored. It induces a shift in the origin 
of 2; = y — a; by an amount /31n(M*/M), which is less than unity since /31n(M*/M) is expected to be of 0(1) or 
smaller. In many cases the local property of fes{z, M) is not significantly changed by this shift. 
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being unaffected. The physics behind this result is both simple and general: the 1/g factor results 
from a probability force in the M direction making the most probable universes in the observer 
region those close to Mc max; while the factor of (3 arises from the shape of the observer boundary 
near Mc, max, as given in Eq. (105). 

We conclude that landscapes having a critical value of M, above which electroweak symmetry 
breaking is impossible, very easily lead to a little hierarchy with z suppressed by both a loop 
factor, (3, and a probability force factor, for example 1/q. Unlike the case of sections 8.2 and 8.3, 
large q does not lead to a large hierarchy. Instead it acts to make M close to M^^max and, since 
^max(^c,max) = 0, this Icads to a little hierarchy. Note that the loop factor can provide an extra 
hierarchy corresponding to 1 - 2 orders of magnitude of fine-tuning, on top of the 1/q factor. 
Hence, this mechanism typically yields M in the TeV region or larger. 

If the physics at M is supersymmetric, this mechanism is the one considered in Ref . [8] . If the 
only scanning parameters are a uniform scaling of the overall scale of supersymmetry breaking and 
of the supersymmetric Higgs mass parameter the scale of Mc corresponds to the scale where 
the determinant of the Higgs boson mass matrix, det(A^|^), passes through zero. To avoid a large 
hierarchy, some other parameter, such as Aqcd should scan. One again needs to assume that, 
for the landscape as a whole, there is some maximum value of the scale where det(Al|^) passes 
through zero, and this is taken to be Mc,max- In this subsection, we have stressed the generality of 
the idea, and that the loop suppression factor appearing in Eq. (107) does not depend on whether 
the physics at M is supersymmetric or not. In the supersymmetric case, fi is determined by the 
specific form of the renormalization group equations for the Higgs boson mass parameters, and 
one finds that j3 ~ 0(0. 1). Nevertheless this is an important accomplishment, since without the 
landscape the ratio between v and M involves a large logarithm ln(M*/f) (in the case of high 
scale supersymmetry breaking) so that A ^ f3\n{M^/v) ^ 0(1), giving 

where Mz is the Z boson mass. The factor of j3, together with 1/q, can easily make M at the 
TeV scale or larger. 

8.6 Summary and discussion 

The conventional naturalness argument implies that the mass scale M of new physics beyond the 
Standard Model that generates electroweak symmetry breaking should not be much larger than 
the Higgs vacuum expectation value v. In contrast, we have shown that environmental selection 
for nuclear stability can lead to M substantially larger than v. We have considered a very general 
framework that is independent of the model of the new physics, assuming only that it generates 
both positive and negative contributions to the Higgs mass-squared parameter, m\ = {x—y)AM'^. 
If the landscape docs not possess a maximum critical mass Mc^max significantly less than M*, the 
relevant observer boundaries are shown in Fig. 17 when only x, y and M scan, and in Fig. 18 
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when yu,d,e, oc and Aqcd also scan. In both cases we find 

• A large hierarchy is generated by a distribution for M (M) that grows sufficient with M 
(M) to overcome the narrowing observer region of Fig. 17 (Fig. 18). While this does not 
require a very strong peaking, the distribution must continue growing over many orders of 
magnitude in M (M). 

• A little hierarchy is generated when the distribution for M (M) does not grow sufficiently 
to generate a large hierarchy, and the distribution for z is strongly peaked to low values of 
z. This strong peaking need only persist for one or two orders of magnitude in z. 

In both cases there is necessarily a large cut factor, either from a mild growth in the distribution 
for M (M) over many decades for a large hierarchy, or from a strong distribution for z over a 
much more limited range for a little hierarchy. In general, a significant variation of a distribution 
function can arise from many sources — the distribution of vacua in the landscape, the population 
of these vacua, integrating out parameters, and the number density of observers — so that in the 
multiverse the existence of a hierarchy, either little or large, is not surprising. 
For landscapes with a maximum critical mass M^max 

• A little hierarchy is generated by a distribution that favors large M near Mc,max- The size 
of the little hierarchy depends on both the strength of this distribution, and also on a loop 
factor that arises from the beta function ior x — y. The combination of these factors makes 
it probable that M is as large as several TeV. 

In section 10 we argue that, in the presence of WIMP dark matter, a strong distribution for M 
can arise from the probability distribution for the cosmological constant. 

Until now we have not specified the physics behind x and y; here we consider a few simple 
schemes, stressing how a cancellation between x and y could lead to a strong distribution for 
z and hence a httle hierarchy. The quadratic divergences of the Standard Model contribute to 
both X (e.g. the SU{2) gauge contribution) and y (e.g. the top quark contribution). At mass 
scale M suppose that these quadratic divergences are cut off by particles of mass Mw' and Mf , 
respectively. If the theory at M is supersymmetric, these are the wino and top squark masses, 
while if the theory at M is non- supersymmetric (for example, with composite Higgs dynamics) 
they are some states of the model. In all of these schemes, integrating out the physics at M gives 
a low energy theory with 

< ~ r£,M^. - ^M,l (108) 

where g is the Standard Model SU (2) gauge coupling. The numerical coefficients are model depen- 
dent and extra contributions to are expected, for example those that regulate the hypercharge 
and Higgs-quartic divergences; but neither of these affects the arguments below. 
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Whether the theory at M is supersymmetric or not, it may well be that M arises as a dimen- 
sional transmutation, in which case it is reasonable that the distribution for In M is sufficiently 
fiat that a large hierarchy does not develop. In this case, what distribution for the parameters in 
Eq. (108) would lead to a little hierarchy? 

We need a distribution fz{z) peaked at low values, i.e. in Fig. 17 or Fig. 18 the probability force 
must have a large component downwards. Since z = is not a special point from the fundamental 
theory point of view (sec the discussion in section 2.2), this implies that most universes will be 
in the phase with z < 0, i.e. m| > 0. In the multiverse, the positive term in Eq. (108) typically 
dominates over the negative term. Such a distribution could arise in several ways. For example, 
suppose that Mw/Mf does not scan, but g and yt do. If g has a distribution peaked at a (much) 
larger value than yt, then most universes will have negative z, and near the observer boundary 
the probabihty force will be towards smaller values of z. The shapes of the distributions for g 
and yt need not be power law, as assumed for simplicity in section 8.4. For example, they could 
be Gaussians with peaks at g and yt, with g sufficiently larger than yt, so that most universes 
have > 0. The few universes that have m| negative will typically have low z and therefore 
a little hierarchy. For narrow Gaussians for g^ and yf with the standard deviations Sg and 5y, 
respectively, we find 



where Zc = {Ag/Ay)g'^ ~ > and 6"^ = {Ag/ Ay)'^6g + 6^, with Ag^y positive coefficients defined 
by ml = {Agg'^ — Ayy^)M'^. Alternatively, it could be that the scanning of the masses Mw/ and 
Mf is more important than that of the couplings, and that the little hierarchy results because 
the distributions typically give Mw' (much) larger than Mf. 

In particular models, it is possible to see other situations that lead to a little hierarchy. For 
example, in the minimal supersymmetric standard model it could be that the distributions for 
the supersymmetric Higgs mass parameter /i and the scale of the soft supersymmetry breaking 
mass parameters m differ. If /i is typically (much) larger than m, then most universes do not 
have electroweak symmetry broken by (h) ^ 0. This generically leads to a strong distribution 
preferring low z, and universes in the observer region will have a little hierarchy. More generally, 
in the minimal supersymmetric standard model the electroweak phase boundary takes the form 



with the soft Higgs mass-squared parameters m|^^ and m^^^ depending on other parameters of 
the theory via renormalization group scaling. Any multiverse distribution for the parameters of 
the model that typically makes the left-hand-side of Eq. (110) larger than the right-hand-side 
will generically lead to a little hierarchy. The scale of the superparticle masses are then raised 




(109) 



(i/^r+<j(i/^r+<j = i/^i?r, 



(110) 
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significantly above v and the measured values of the parameters should be close to satisfying the 
critical condition, Eq. (110). 



9 Electroweak Symmetry Breaking as Observer Boundary 

In the last section we assumed that the relevant observer boundaries for selecting the electroweak 
vacuum expectation value, (h) = v, were those of neutron, deuteron and complex nuclei stability. 
Of these three boundaries, the requirement that some complex nuclei arc stable seems clearly 
to be the most robust requirement for observers. For example, if the neutron is stable, nuclear 
energy is produced in diffuse protogalaxies rather than in stars. While this is a drastic change 
from our universe, some form of observers might still be possible. On the other hand, it is harder 
to imagine that some complex observers develop in the world in which the only stable nucleus is p 
or In this section, we retain only the complex stable nuclei boundary, dropping the neutron 

and deuteron (in) stability requirements from the observer boundary. This allows much smaller 
values for v. How small can v become while remaining in the observer region? If electroweak 
symmetry is broken dominantly by the QCD condensate, there is a strong washout of the baryon 
asymmetry of the universe due to sphaleron effects. This implies that complex structures involving 
baryons do not arise ii v ^ ^qcb, so that there is a complexity boundary for v near Aqcd which 
is therefore a candidate for being part of the observer boundary. Since Aqcd is much smaller 
than the scale of the electroweak symmetry breaking sector, M, we can speak of this boundary 
as the phase boundary between {h) = and (h) ^ phases. In this section, we use only two 
boundaries: complex nuclear stability and the Higgs breaking of electroweak symmetry, which 
were both classified as catastrophic boundaries in section 4. 

We consider the generic electroweak symmetry breaking sector discussed in the previous sec- 
tion, with relevant scanning parameters x, y and M and a hierarchy 



With yu,d,e, Oi and Aqcd fixed, the relevant observer region is now as shown in Fig. 22. Compared 
to Fig. 17, the observer boundary at f+ remains while that at V- has been replace by one at 
V ~ Aqcd- This greatly enlarges the observer region, with v varying over about 4 orders of 
magnitude compared to a factor 4 in section 8. Despite this, we find that there is very little 
change in the physical picture of environmental selection. 

For example, in the case that y is fixed and x and M have polynomial distributions as 
in Eq. (90), the effective distributions for z and v are not changed and given as before (see 



Eqs. (92, 98, 99)) by 




(112) 
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Figure 22: Sketch of the 2- dimensional observer region in the M-z plane. The two boundaries of 
the observer region depicted are v = v+ and v ~ Aqcd- 



Once again, q > 2 gives a large hierarchy and ^ < 2 a little (or no) hierarchy. Furthermore the 
size of the httle hierarchy is largely governed by n, as discussed in section 8.4. One difference 
is that the range of v in the observer region is now much larger; how reasonable is it that our 
universe is about a factor of 5 from v+ and 3 orders of magnitude from Aqcd? For ^ q ^ 1 
the observed value of v is quite typical, but for larger q, v becomes progressively more peaked 
near v+. However, this cannot be viewed as evidence against a large hierarchy. For large q it is 
unlikely that Eq. (113) is the correct distribution near as v approaches this observer boundary 
so successively more nuclei become unstable, so that the observer factor, n, is likely to become 
smaller. 

Given the observer region of Fig. 22 we can understand the origin of large or little hierarchies 
from the qualitative features of the distribution f{z,M), without recourse to any particular 
functional form. The observer boundary at f+ causes a narrowing of the observer region at large 
M, so that a large hierarchy results only if the distribution gives a sufficient preference to large 
M to overcome this narrowing. If the probability force in the M direction is insufficient, then 
the size of the little hierarchy depends on the strength of the probability force in the negative 
z direction. At first sight a force in the negative z direction apparently pushes the electroweak 
vacuum expectation value close to the value Aqcd at the baryon washout boundary. This is 



and 




for q>2, 
for q <2. 



(113) 
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incorrect: such a force makes the most probable region that with low but the value of v that 
this corresponds to changes with M. This is why Eq. (113) involves q and not n. 

We again conclude that the little hierarchy problem (supersymmetric or not) is very easily 
solved by environmental selection. Strongly varying distribution functions of the electrowcak 
symmetry breaking sector are able to generate either a large or little hierarchy; furthermore this 
is not sensitive to which observer boundary is used to limit the lower value of v. 

10 Connections to the Cosmological Constant 

In this paper we argue that evidence for the multiverse can be found in unnaturalness in the 
cosmological constant, nuclear physics, and electroweak symmetry breaking. This implies that 
scanning of parameters is occurring in all three different arenas. Is it correct to consider each 
scanning problem separately, or should there be a combined treatment? In section 8.3 we provided 
an analysis of the combined scanning in the Standard Model and electroweak symmetry breaking 
sectors. In this section we study connections with the scanning of the cosmological constant. 

In section 6.2, we argued that "integrating out" a set of parameters Xf, modifies the distribution 
function for a set Xa 

fes{Xa) = f priorix a, Xb)dXb, (114) 

if the observer region O for Xb depends on Xa- In the present case we consider Xb = A. and Xa to 
be the set of scanning parameters of the Standard Model and the electroweak symmetry breaking 
sector. Assuming the multiverse probabihty distribution for A to be flat on a linear scale 

— /prior (^^a) PNhiXa): (115) 

where is the energy density of the imiverse when it goes non-linear. Hence the distribution 
function for the parameters of the Standard Model and the electroweak symmetry breaking sector 
is modified by the scanning of the cosmological constant if pnl depends on the parameters Xa- 
In our previous analyses, the assumed form of the distribution should apply to fcs, rather than 
to /prior- This implies that the scanning of the cosmological constant can affect the form of the 
distribution functions that appeared in the analyses in previous sections. 

In this section, we explore some consequences of a nontrivial dependence of pnl on Xa- Perhaps 
the simplest origin for such a dependence is the case of WIMP dark matter. The temperature 
of matter-radiation equality, Tgq, depends on the WIMP mass via its annihilation cross section, 
Teq oc 1/(Ja ~ "^wiMP- Hcncc Pnl ~ Q'^^eq ~ Q^^wiMP) where Q is the primordial density 
perturbation. Working with the generic electroweak symmetry breaking sector introduced in 
section 8, with mass scale M and dimensionless parameters Xi, the WIMP mass is proportional 
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to M so that 

Teq = g{xi) (116) 

where g{xi) > is a model dependent function, depending on the mass and interactions of the 
WIMP. The effective distribution function for the electroweak symmetry breaking sector now 
becomes 

/PNL(a:i,M) 
dA ~ fpriorixi, M) g{xiYM^Q\ (117) 

Note that here we have assumed that the parameter Q does not scan. If Q also scans, the result 
of integrating out cosmological parameters is altered from Eq. (117), as we will see below. 

If Q does not scan, Eq. (117) shows that the effective distribution for the parameters of 
the electroweak symmetry breaking sector, {xj, M}, can acquire an important component from 
integrating out the cosmological constant. In particular, the cosmological constant may provide 
a strong probability force to larger values of M through the factor. At first sight it appears 
that this factor drives a large hierarchy, M ^ v, but this is not the case. A crucial issue is: what 
stops the runaway to large Ml It is important to separate two cases. 

In the first case, the factor ~ {g[xi)M'^}'^ in /efj pushes the amount of dark matter up 
to the maximum allowed by some astrophysical observer boundary, for example that of stellar 
collisions in galaxies [33], so that Teq is essentially fixed to the boundary value Teq = T^f^,*- The 
distribution function for the parameters relevant for electroweak symmetry breaking, {x, A/i, M} 
where x, Xh C Xi, is obtained by integrating out parameters that appear in g{xi) but not in m| 
or Xh, within the observer region Teq < Teq,*. This in general leads to a complicated dependence 
of fcs on {x,y, Xh, M}, not in the simple form of Eq. (117). Note that as M grows beyond the 
TeV scale, since Tgq stays to be Teq,*, cancellations must occur in g{xi) so that mwiMP ^ M. 
This implies that a large hierarchy can develop only if /prior has a strong preference towards larger 
values of M. 

In the second case, the most probable region of scanning parameter space does not lead to 
the dark matter density being on the edge of its maximal value determined by astrophysics. The 
values of the parameters {xi, M} are determined by other physics, including electroweak symmetry 
breaking. The distribution of the parameters relevant for electroweak symmetry breaking may 
take a form f^s ~ fT^q, so that the factor of g{xiYM^ can play an important role in electroweak 
symmetry breaking. Some examples of this are given in the following subsection. 

A completely different situation arises if the space of scanning parameters is increased to 
include Q. As illustrated in Fig. 23, the probability force from the cosmological constant no 
longer acts in the direction of increasing mwiMP, but in the direction of increasing pnl, so that 
now the issue becomes: what stops the runaway in the pnl direction? It is again an astrophysical 
limit on the amount of dark matter, but this is now a limit on pnl rather than on Teq: pnl < Pnl,*- 
The probability distribution along this astrophysical boundary now depends on the distribution 
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Figure 23: The probability force from the cosmological constant in the Q-Tcq plane. The contours 
of pnl are drawn by dashed lines, while the observer boundary from astrophysics, pnl = Pnl,* is 
depicted in the solid line. 

fgiQ). Integrating first A and then Q over the observer region, we obtain 

U{x^, M) OC /pHor(x„ M) j^''''^''' '^"''^ ' Jq^Q) dQ , (118) 

where Tgq is given by Eq. (116). The induced contribution to the effective distribution now depends 
on the form of the original distribution for the density perturbation Q, since the observer boundary 
for Q depends on Teq = T^q^Xi, M). However, in the simple case that /g oc 1/Q, so that the Q 
distribution is flat on a logarithmic scale, fes{xi,M) ~ /prior (a^i, Af), and the TJ^ contribution 
is removed. Hence, a nontrivial contribution to f^s is generated by a nontrivial distribution for 
InQ, not by the probability force from the cosmological constant. 

10.1 Examples of the cosmological constant affecting electroweak sym- 
metry breaking 

We have seen that integrating out cosmological parameters with WIMP dark matter can gener- 
ate a nontrivial component in the distribution of parameters relevant for electroweak symmetry 
breaking. In general the induced distribution takes a complicated form depending on what pa- 
rameters scan and how they are determined. In the case that Q does not scan and the dark matter 
density is not determined by an astrophysical bound, the effective distribution function receives 
a factor of g{xi)'^M^, as shown in Eq. (117). Assuming that parameters appearing in g{xi) but 
not in ml or are determined independently of M, the distribution for the parameters x, y and 
M, appearing in m\ = {x — y)AM, takes the form 

/efr(a;,y,M) ~/pHor(a;,2/,M)^(x,|/)V, (119) 
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where g{x,y) is related to g{xi) in Eq. (116). In this subsection we illustrate how the scan- 
ning of cosmological parameters can affect electroweak symmetry breaking, using the example of 
Eq. (119). 

Suppose that /prior is a rather mild function of x, y and M. The effective distribution fes is 
then determined essentially by the factor g{x,y)^M^. The fact that the dark matter density does 
not saturate the astrophysical bound implies that the runaway to larger M due to the factor 
should be stopped not by the condition T^q < Tcq,* but by some other physics. The simplest 
possibility is that M does not scan in the multiverse. At the end of section 8.2 we argued that 
if M is fixed then a large hierarchy is to be expected, as in Eq. (82). However, this conclusion is 
reversed if Aqcd scans. Since the nuclear physics boundaries depend on M only through the ratio 
M/Aqcd, it makes no difference to the argument on the size of the hierarchy whether M scans 
or Aqcd scans. Allowing x, y and Aqcd to scan, the observer region is that of Fig. 17, with M 
replaced by 1/Aqcd- Taking a mild distribution for Aqcd avoids a large hierarchy. ^° The amount 
of the little hierarchy is then determined by the factor g{x, y)^. For example, for g r-> x"^ and a 
logarithmic distribution for Aqcd, (z) ~ l/(4m-|-2), so that a little hierarchy of order 10~^ - lO"-*^ 
can be easily obtained for m a factor of a few. In general, with gi^ strongly preferring x > y, 
we obtain a little hierarchy. The crucial physics here is that the large scale structure observer 
boundary, A ^ pnl, depends sensitively on the electroweak symmetry breaking parameters x and 
y through WIMP dark matter. The probability distribution for the cosmological constant then 
generates an effective distribution for x and y, which can be sharply varying. 

Another possibility of preventing the runaway to large M arises if the range for the scanning 
of X and y are such that there is a maximum energy, Mc^max, above which the ml parameter is 
always positive. As shown in detail in section 8.5, a strong probability force to large M leads to 
a little hierarchy, providing the nuclear observer boundaries of v± are able to scan, for example 
via a scanning of Aqcd- Here wc simply point out that the strong force in the M direction can 
arise from integrating out the cosmological constant, i.e. from the factor in Eq. (119). This 
leads to a little hierarchy without the need for any cut factor beyond that for the cosmological 
constant. In the case of supersymmetry, with dj\f ~ M^dM {q a; 0(10)) the multiverse allows an 
improvement in the naturalness by a factor oi ^ q/ P !v O{100). 

11 Conclusions 

Environmental selection on a multiverse is a radical departure from conventional methods of 
fundamental physics for explaining physical phenomena. Nevertheless, in light of the cosmological 

^°In fact, consistency of the setup requires that the hierarchy is not very large, since a value of M much 
larger than the weak scale makes the natural size of Teq ~ SttM^/Mpi much larger than Teq,*, implying that the 
astrophysical bound is saturated. 
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constant problem and the discovery of dark energy, it warrants further exploration. Will sufficient 
evidence emerge to convince us that the multiverse exists? 

Since we cannot directly explore other universes, it may be questioned whether evidence for 
the multiverse can be found at all. However, theories in physics and cosmology that cannot be 
directly tested in the laboratory are far from new — one only has to think of unified theories and 
infiation. The inability to make direct laboratory tests of the new particles and interactions does 
not put these theories beyond the realm of science; rather, it leads to careful investigations of 
whether they make successful indirect numerical predictions for data that cannot be satisfactorily 
explained by other means. Thus for any new theory or framework that cannot be directly probed, 
two questions are important 

(a) In the absence of the new theory, is there numerical data that cannot be adequately explained 
using other known theories? Indeed, does the data represent a problem for existing theories? 

(b) Does the new theory provide a numerical understanding of the data, thus solving the previous 
problem? 

For inflation, the data of (a) includes the flatness and isotropy of the universe, and the spec- 
trum of density perturbations. So far, these highly signiflcant cosmological problems have only 
been solved by infiation. Similarly, unified theories provide an understanding of gauge coupling 
constant unification. If the numerical understanding of the data is sufficiently precise, and thought 
to be better than competing theories, then the new theory may become the provisional standard 
view. As long as the numerical significance of its predictions is not overwhelming, one must stress 
the provisional nature of the understanding and the importance of seeking both further arenas 
in which it can be tested and new competing theories. For example, proton decay would provide 
further evidence for unified theories, but despite intensive searches, such evidence is still lacking. 

Seeking evidence for the multiverse is no different in principle than seeking evidence for other 
theories that cannot be directly probed in the laboratory. In this paper we have argued that 
evidence for the multiverse can be found in three different arenas: the cosmological constant, 
nuclear physics and electroweak symmetry breaking. In all three cases the conventional approach 
based on symmetries has not provided a numerical understanding of the data, rather in each case 
it leads to naturalness problems. The observed values of parameters are very close to special values 
that are critical for the formation of some complex structure, yet this closeness is not adequately 
explained by any symmetry. An observer region in the parameter space of some theory is defined 
by requiring the existence of certain complex structures necessary for observers. Unnaturalness 
results if the observer region is very small compared to the entire parameter space, or if we observe 
values of the parameters very close to the boundary of the observer region. We have introduced a 
naturalness probability, P, that allows a numerical evaluation of these problems. The value of P 
can be highly dependent on the theory T under consideration. The evidence for unnaturalness. 
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(a), is governed by the maximal value of Pt that can be obtained in simple theories: the lower 
the maximal Pt, the more severe the naturalness problem. 

In each arena the multiverse easily and generically solves the naturalness problem. Each 
arena is somewhat different, and we summarize our results for each below; but there are also 
some common features. In all three cases, environmental selection elegantly explains why our 
universe is not to be found in the largest region of parameter space. Furthermore, the multiverse 
distribution is likely to have a strong dependence on the parameters of the low energy effective 
theory, through the landscape of vacua of the fundamental theory, the population mechanism, 
integrating out parameters, and from the physics that determines the density of observers. With 
a strongly varying distribution, it is most probable to observe a universe close to the observer 
boundary, solving naturalness problems and leading to predictions. Given the current theoretical 
status, the multiverse appears to us to provide the most elegant and plausible prediction for 
several parameters, including A and mu,d,e- 

The cosmological constant problem is the most severe naturalness problem, with a naturalness 
probability in the range 

Pa ^ (10-^2° - 10-^°). (120) 

The robustness of the observer boundary is particularly convincing — how could observers form 
in a dilute gas of inflating elementary particles? The interesting open questions are whether 
our universe is sufficiently close to the observer boundary, and how runaway behavior can be 
prevented if Teq and Q scan. These questions, however, are secondary: a notoriously intractable 
problem has an elegant solution that predicts dark energy. 

Astrophysicists have often remarked on the special values of parameters required for a variety 
of phenomena in nuclear physics. We have obtained the observer region in the 4-dimensional 
parameter space mu,d,e/ ^qcDi ol resulting from the stability boundaries for neutrons, deuterons 
and complex nuclei, as shown in Fig. 5. We find this observer region to be small, with our 
universe within (10 - 30)% of the neutron stability boundary, so that in the Standard Model 
the naturalness probability is Pnuc.SM ~ (10^^^ - lO^''). Even if we could construct a theory 
of flavor with successful, precise predictions for the Yukawa couplings yu4,e and at the same 
time find a theory that correctly predicted the weak scale, there would still be a naturalness 
probability associated with the value of Aqcd that we estimate to be ^ (0.01 - 0.04). Despite 
decades of experiments on fiavor physics, with recent increasing levels of accuracy in B meson and 
neutrino physics, progress on a theory of fiavor has been limited. The most promising theories 
appear to be based on fiavor symmetries, with a sequential pattern of transferring symmetry 
breaking to successive generations of quarks and charged leptons. There are a great number of 
candidate theories, but none is sufficiently promising to be widely recognized as the standard. 
We have estimated that such a lack of progress in fiavor physics, especially in the first generation 
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masses, decreases the naturalness probability of the nuclear observer region by about an order of 
magnitude, leading to 

Pnuc^(10-'-10-')- (121) 

We think that this estimate is conservative. For example, in theories with Abelian flavor symme- 
tries this estimate does not take into account that there are many simple ways of assigning charges 
to the three generations. Given the state of theories of flavor, we suspect that the naturalness 
probability of the nuclear observer boundary is less than this conservative estimate. Suppose that 
mu,(i,e bad each been significantly different, for example by an order of magnitude. It would still 
be possible to accommodate this in theories with Abelian flavor symmetries, with about the same 
level of success as for the actual observed values, by changing the charges of the first generation. 

The nuclear naturalness problem certainly requires more than 1% fine-tuning, and has not 
received sufficient recognition. If symmetries rule flavor physics, it must be viewed as purely 
accidental; while in the multiverse it can be viewed as a prediction. After decades of studying 
theories of flavor, the current status is that we arc far from a convincing explanation for the masses 
of the electron, the up quark and the down quark. The multiverse allows simple arguments that 
relate these masses to the QCD scale. For example, if the multiverse distribution strongly favors 
isospin restoration, then rrie and rrid — m„ are expected to be close to 5em — 1.0 ± 0.5 MeV, 
the electromagnetic mass difference of the proton and neutron. It is true that such multiverse 
predictions require an assumption on the form of the distribution function; but such assumptions 
may be much simpler than the choice of flavor group, representations and sequential symmetry 
breaking of the standard approach. 

The multiverse predictions for mu,d,e/ ^q,CT> can be retained even if unified and/or flavor sym- 
metries describe the overall pattern of quark and lepton masses and mixings. This allows us to 
preserve particular successful relations, such as ~ sffn^fnis or m\, — rrir at the unifled scale. 
The only requirement is that the theory contain three independent scanning parameters that 
allow mn,d,e/AQCD to be environmentally selected. 

A simple, natural theory of electroweak symmetry breaking is lacking. The natural regions of 
simple technicolor, supersymmctric and composite Higgs models have been excluded. Precision 
measurements of electroweak observables, together with direct searches for the Higgs boson, have 
led to successive increases in unnaturalness, leading to 

Pbwsb ^ (10-' - 10-^), (122) 

in simple models. This problem in the data is of type (a) and is widely appreciated, so that much 
recent research has focused on building models that alleviate the problem. We have shown that, 
no matter what the ultimate physics at mass scale M behind electroweak symmetry breaking, 
environmental selection on a multiverse leads very easily to M ^ v, with either a little or 
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large hierarchy. This is a very robust result requiring only that some parameters of electroweak 
symmetry breaking scan, and that the multiverse distribution is strongly varying. A distribution 
/(M) ~ M*^ gives a large hierarchy if g > 2 and a little hierarchy (or no hierarchy) for q < 2, 
which includes the important case of M being induced by a dimensional transmutation (g = 0). 
For q < 2 the little hierarchy, f^/M^, gains a factor of ^ 1/n from the multiverse, where n 
describes the peaking of a distribution in some other parameter of the electroweak symmetry 
breaking sector. In multiverses where electroweak symmetry breaking is only possible if M is 
below some critical value, Mc^max, a little hierarchy develops from a distribution favoring large M. 
The size of the hierarchy, w^/M^, is enhanced by a loop factor j3, gaining a factor of ^ so that 
M is typically in the TeV region or larger. These results are independent of what other Standard 
Model parameters are scanning, including Yukawa couplings and Aqcd, and do not even depend 
on whether selection is happening at the nuclear stability boundaries, or at the phase boundary 
for electroweak symmetry breaking itself. Hence we stress that the multiverse provides a very 
general solution to the hierarchy problem, whether little or large. 

All three naturalness problems have the common feature of being solved by a multiverse distri- 
bution that makes observers typically close to an observer boundary. There may be connections 
between the three problems, arising from integrating out certain parameters in a more funda- 
mental theory. For example, the probability force driving a little hierarchy may originates from 
the probability distribution for the cosmological constant, with WIMP dark matter acting as a 
mediator. 

The current status of naturalness in electroweak symmetry breaking, Eq. (122), indicates that 
either we have not yet arrived at the right theory or that environmental selection is playing an 
important role. The LHC will determine the correct interpretation of Eq. (122), leading us either 
to a new natural theory, or to a third arena for multiverse evidence. Even if the naturalness 
probability is much larger than for the cosmological constant, the pervasive pattern of a finely 
tuned universe will make the multiverse much harder to dismiss. The importance of the LHC 
in this regard cannot be overemphasized: for the cosmological constant and nuclear naturalness 
problems we are stuck — we may not be able to experimentally determine the relevant theory 
beyond what we already know, in which case increasing the naturalness probability is a theoretical 
enterprise. However, the LHC will teach us a great deal about the theory of electroweak symmetry 
breaking and hence will lead to a better determination of Pewsb- One possibility is that the LHC 
will reveal a completely natural theory that we have not been able to invent. Below we mention 
a few examples where LHC data could determine a small value for Pewsb- 

If the LHC discovers a light Higgs boson of mass mniggs and sets a hmit of Mcoi on new colored 
particles, then 




(123) 
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This result is true in the vast majority of theories, although some counterexamples are known. 
Another possibility is that only a light Higgs boson is discovered, with a mass very close to 
the vacuum instability limit of the Standard Model. While Eq. (123) still applies as a direct 
consequence, there is the additional implication that the hierarchy is large with a very much 
smaller Pewsb- Evidence for a large hierarchy could also emerge from the discovery of split 
super symmetry. Alternatively, the discovery of weak scale supersymmctry, with a light Higgs 
boson and a top squark heavier than ^ 1 TeV would indicate a little hierarchy with a naturalness 
probability 



where Mmess is the messenger scale of supcrsymmetry breaking. 

The discovery of dark energy has verified a remarkable prediction of the multiverse; but this 
could be undermined by the discovery of an alternative solution to the cosmological constant 
problem. The nuclear stability boundaries imply at least 1% fine-tuning in any known theory, 
and the multiverse allows a striking understanding of mu^d,e- Our current theories of electroweak 
symmetry breaking are unnatural; a confirmation by the LHC would solidify evidence for the 
multiverse in a third arena. Even with Pnuc and Pewsb much larger than Pa, the three arenas 
together would provide significant, robust evidence for a multiverse. 
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